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RELATIVISTIC  CONTINUUM  MECHANICS 
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ABSTRACT 

The  basic  lave  of  a  special  relativistic  theory  of  continuous 
media  suitable  for  the  treatment  of  electromagnetic  interactions 
with  materials  are  formulated.  The  kinematics,  dynamics  and  thermo¬ 
dynamics  of  a  continuum  are  dis~ussed  from  a  relativistic  viewpoint. 
Constitutive  equations  are  deduced  for  thermoeiastic  solids,  thermo- 
viscous  fluids  and  electromagnetic  materials. 
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INTRODUCTION 


In  this  article  there  is  presented  a  nonlinear  relativistically 
(special)  invariant  theory  of  continuous  media.  The  object  of  such 
a  study  is  a  consistent  treatment  of  the  interaction  of  electromag¬ 
netic  fields  with  the  deformation  of  matter.  In  recent  years  several 
theoretical  works  have  appeared  concerning  the  simultaneous  action 
of  large  deformations  and  electromagnetic  fields  cn  material  bodies. 

In  general  these  researches  have  either  considered  only  static  de¬ 
formations  (cf.  Toupin  [1],  Eringen  [2],  Jordan  and  Eringen  [3],  [U]) 
or  developed  a  dynamical  theory  along  nonrelativistic  lines  (cf. 

Toupin  [5],  and  Dixon  and  Eringen  [6]). 

It  is  well  known  that  the  invariance  group  of  the  basic  equations 
of  electromagnetism  is  the  Lorentz  group.  The  modern  theories  of 
continuum  mechanics,  however,  make  use  of  the  Galilean  group  for  the 
basic  laws  of  motion  and  the  invariance  under  the  group  of  rigid 
motions  for  the  constitutive  theory.  At  the  turn  of  the  century,  the 
invariance  of  the  laws  of  mechanics  under  the  Galilean  group  was  dis¬ 
carded  at  least  for  the  physical  phenomena  which  fall  within  the  scope 
of  relativistic  considerations.  Thus  we  believe  that  a  satisfactory 
and  consistent,  theory  of  electromagnetic  interactions  with  deforming 
materials  cannot  be  obtained  until  the  ground  rule  (the  invariance 
principles)  for  mechanics  and  electromagnetism  is  taken  to  be  the  3ame. 
Within  the  scope  of  the  special  theory  of  relativity  the  most  natural 


and  the  simplest  invariance  principle  is  the  Lorentz  group  of  trans- 
formations . 

Following  the  scheme  of  classical  theories  of  continuum  mechanics 
Eringen  [7],  Truesdell  and  Toupin  [8],  we  begin  our  study  with  the 
geometrical  and  kinematical  description  of  deformation  and  motion. 
Afterwards  the  basic  lava --concerning  the  motion  and  physical  phenomena 

are  introduced.  The  formulation  is  then  completed  with  the  const! - 

i 

tutive  theories  in  accordance  with  certain  methodological  principles 
set  forth  in  [9]. 

In  the  literature  the  kinematics  of  continuous  media  is  seldom 
discussed  from  a  relativistic  viewpoint.  Exceptions  are  the  works  of 
Bressan  [lO]  and  Toupin  [llj.  BreBsar.  has  formulated  a  relativist! - 
cally  (general)  invariant  kinematics  by  noting  that  the  motion  of  a 
material  particle  with  undeformed  coordinates  X  (K  =  1,  ...,  3)  can 
be  described  by  x*1  .  x*1  (X*  t)  (n  =  1,  4)  where  t  is  some 

temporal  parameter.  Observing  that  such  a  temporal  parameterization 
is  highly  arbitrary,  he  stipulates  that  the  description  of  the  de¬ 
formation  is  independent  of  the  choice  of  t  . 

(Xir  approach  to  kinematics  is  based  on  an  extension  of  the  work 
of  Toupin  [ll].  ( Xir  strain  measures  are  invariant  under  the  Lorentz 

group  and  they  coincide  with  thoss  of  Bressan  under  appropriate  modifi 
cations.  The  world  velocity  that  is  selected  is  a  kinematical  one. 

In  several  works  the  world  velocity  of  a  continuous  body  is  defined  as 
either  the  time -like  eigenvector  of  the  energy -momentum  tensor  (cf. 
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Synge  [12]  or  Lichnerovici;  [13])  or  as  the  unit  vector  parallel  to  the 
momentum  density  (cf.  Miller  [l4]  or  Landau  and  Lifehitz  [ 15 3 • ■ 

Miller  actually  employs  also  the  kinematical  velocity. )  This  velocity 
of  energy  propagation  will  coincide  with  our  kinematical  velocity  only 
in  special  cases,  for  example,  when  the  hc»t  conduction  and  electro¬ 
magnetic  phenomena  are  neglected. 

For  basic  mechanical  and  thermomechanical  balance  conditions 
ve  eelect,  in  Chapter  II,  the  conservation  of  particle -number,  the 
balance  of  energy -momentum,  the  balance  of  moment  of  energy -momentum 
and  the  second  lav  cf  thermodynamics.  We  consider  the  conservation  of 
particle  number  as  the  appropriate  generalization  of  the  conservation 
of  mass  in  nonrelativistic  continuum  mechanics.  Thie  point  of  view 
ie  accepted  by  several  authors  (for  example,  Landau  and  Lifehitz  [15], 

Van  Lantzig  [l6]  and  Eckart  [ 1 Y ] ) •  An  alternate  viewpoint  ie  that  the 
classical  conservation  of  mase  is  obtained  through  the  equation  of 
energy  in  the  limit  as  the  epeed  of  light  approaches  infinity  (for 
exantpie,  Bergmann  [l8],  Thomas  [19]  or  Edelen  [20]).  From  thie  point 
of  view  it  ie  difficult  to  consider  thermodynamics  without  an  additional 
axiom.  In  nonrelativietic  theories  the  thermodynamical  balance  equations 
are  deduced  from  the  energy  equation.  It  would  seem  appropriate  to 
do  the  earns  in  relativistic  theories.  Thus,  for  example,  the  energy 
equation  should  be  ueed  to  derive  the  equation  of  heat  conduction. 

The  form  of  the  balance  of  energy -momentum  ie  well  known  and 
hardly  needs  any  explanation.  It  should  be  noted  that  in  a  relativistic 
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theory  part  of  the  momentum  is  due  to  nonmechanical  sources  such  as 
the  heat  flux  and  stress  tensor.  In  most  of  the  works  on  relativistic 
theories  of  conti  nua  the  energy -momentum  tensor  is  assumed  to  be 
symmetric.  This  is  usually  deduced  from  the  balance  of  angular  mo¬ 
mentum.  In  the  classical  theory  of  continuum  mechanics  the  e trees 
tensor  is  symmetric  only  if  there  are  no  torques  acting  on  the  body 
and  the  body  does  not  possess  an  intrinsic  spin.  Ve  generalize  this 
idea  to  the  relativistic  case  by  introducing  a  spin  tensor  and  a  body 
torque  tensor  in  four  dimensions.  The  idea  of  a  spin  tensor  can  be 
found  in  such  works  as  Bcgoliubov  and  Shirkov  [21]  and  Papapetrou  [22]. 

If  the  spin  and  torque  tensors  vanish,  it  follows  from  the  balance 
of  moment  of  energy -momentum  that  the  energy -momentum  tensor  is  sy;im»etric. 

The  second  law  of  thermodynamics  is  formulated  in  a  manner 
analogous  to  the  law  given  by  Eckart  [IT ]  with  the  provision  that 
it  is  considered  as  a  restriction  on  the  form  of  the  constitutive 
equations  rather  than  on  the  electromechanical  process. 

In  Chapter  III  the  methodological  principles  for  formulating 
constitutive  equations  are  set  down.  They  are  obvious  generalizations 
of  those  employed  in  nonrelativistic  continuum  mechanics.  The  major 
difference  between  these  principles  and  those  of  modern  continuum 
mechanics  is  the  requirement  that  the  constitutive  equations  are  to 
be  form-invariant  under  the  Lorentz  group.  This  is  adequate  for  the 
theories  treated  in  this  paper.  A  different  approach  using  the  idea 
of  non-sentient  response  is  to  be  fc*iM  elsewhere,  Bragg  [25]. 
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As  example b  of  the  above  formulation  of  constitutive  equations, 
relativistic  theories  of  thermo solids  and  thermo viscous  fluids  are 
deduced.  The  second  lav  of  thermodynamics  forces  the  acceptance 
of  the  heat  conduction  lav  of  Eckart  [17]  as  opposed  to  those  of 
Bressan  [2k]  or  Pham  Mau  Quan  [25]  vho  attempt  to  generalize  Fourier's 
lav. 

In  Chapter  IV  the  Minkovski  form  of  the  equations  of  electro¬ 
magnetism  are  formulated  in  the  usual  four  dimensional  form .  In 
Chapter  V  the  interaction  between  electromagnetic  fields  and  matter 
is  treated.  An  Interaction  term  is  vritten  dovn  using  the  analysis 
of  Dixon  and  Eringen  [6].  If  one  wishes  to  treat  the  interaction 
of  electromagnetic  fields  and  matter  one  needs  a  physical  model 
(cf.  Dixon  and  Eringen  [6],  Jordan  and  Eringen  [3],  or  the  appendix 
of  Fano,  Chu  and  Adler  [26]).  The  seemingly  fond  wish  to  use  the 
Minkowski  stress  tensor  or  its  symmetric  part  is  erroneous  unless 
the  material  has  no  polarization  or  magnetization,  in  vhlch  case  this 
tensor  is  symmetric. 

In  Chapter  VI  the  constitutive  theory  of  electromagnetic  materials 
is  set  dovn  for  an  elastic  solid  and  a  viscous  fluid.  Such  effects 
as  heat  conduction,  electrical  conduction,  polarisation,  and  magneti¬ 
zation  sure  included  in  the  theory.  The  consequences  of  the  second  lav 
of  thermodynamics  are  fully  investigated. 

Finally  an  appendix  on  the  invariants  of  tensors  and  vectors  for  the 
Lorentz  group  is  included.  In  special  cases  it  is  shown  that  the  use  of 
results  of  the  three  dimensional  orthogonal  group  is  permissible.  *JMs 
facilitates  greatly  the  construction  of  constitutive  equations  since  a 
reworking  of  the  theory  of  invariants  of  teas or s  and  vectors  for  the 
Lorentz  group  would  be  lengthy  and  tedious. 


Rotation 


In  this  paper  we  use  the  standard  tensor  notation  and  summation 
convention  of  repeated  indices.  The  signature  of  the  Lo recti  metric 
(/**)  is  (+  +  +-),  i.e.  711  «  *  1  t  *  -1  >  and  all 
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0  .  It  is  convenient  to  use  the  following  convection: 


*0076  *  *007^ 

where  e^^  is  the  permutation  symbol  and  define 

G07&  .0075 

<  a  -  e 

3ince  det  (7^V)  *  -1 


C0t6  «fl  Wi  77i  Wi 

1  •  ’  ’  ’  7  'WA 


The  short  hand  notation 


Jo 07 


is  often  employed  to  define  the  "cross  product"  between  two  four 
vectors  a3  ,  ba  ,  e.g.,  ft,  x  is  the  vector 

<*«6>r  • :?®  *„  \ 

By  the  outer  product  notation  ft  &  fe  **  will  nean 


1  aab0 


- ,  -  ■*<*;-  :  mnw.  — - - - - - -  .... .  ,  I,  . .  ■  I  . .  ^ 

I 

Tor  a  second  order  tensor  A^  and.  a  four  vector  ve  Till  Bean  by 
the  notation  &  x  £  the  second  order  tensor 

(ftx^  5  «a67*6*7p 
and  by  £  x  ft  the  second  order  tensor 

(axft)^  =  V*er5‘B 

where  vc  raise  and  lower  Greek  indices  with  the  metric  (7°^)  =  (7  _ ) . 

(X 

If  is  a  second  order  tensor  and  a  a  rector  by  £  ft  we  mean 

the  rector 

For  two  vectors  a°  ,  ba  the  inner  product  ft  •  ft  is  defined  by 

.  a  -u 

a  •  b  *  a  ba 

Units 

For  simplicity  the  speed  of  light  c  is  set  equal  to  unity, 
c  °  1  «  In  the  electromagnetic  section  of  this  article  latlonalized 
natural  units  are  employed.  Accordingly  ve  take  the  dielectric 
constant  of  free  space  «Q  *  1  and  since  c  ■  1  ,  the  permeability 
of  free  space  uQ  •  1  . 
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CHAPttH  I 

RBIATIYISTIC  KM3ttTICS  CF  A  COUTIBUOUB  MEDIUM 

Motion 

Although  soat  account*  on  the  relativistic  treatment  of  the 
dgmud.ee  of  a  continuous  medium  may  be  found  in  various  bvxts  on 
relativity  (cf.  Miller  [lk],  Tolmn  [27],  Bergnann  [18],  Landau  and 
Llfshitz  [15,  28],  Synge  [12] ),  very  little  discussion  has  been 
directed  towards  a  relativistically  invariant  kinematics  of  a  con¬ 
tinuous  body.  Tb»  present  chapter  Is,  therefore,  devoted  to  a  dis¬ 
cussion  of  deformation  and  motion  (kinematics)  from  the  point  of 
view  of  the  special  theory  of  relativity. 

In  classical  continuum  mechanics  the  deformation  of  a  body  is 
described  by  specifying  the  mapping  of  one  configuration  of  a  body, 
B  ,  onto  another  configuration,  Bt  • 


The  configuration  B  ,  called  the  reference  state,  (usually  taken  as 
the  undefonaed  body)  can  be  made  explicit  by  specifying  the  rectangular 
coordinates  X*  ,  (K  •  1,  2,  5)  of  a  set  of  material  points,  (B)  ,  in 
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throughout  body  3  except  possibly  sons  singular  surfaces.  Uses  sad 
points.  In  fact,  unless  otherwise  stated  re  shall  sot  only  as  suae 
the  validity  of  (1.3)  throughout  B  but  also  the  existence  of  partial 
derivatives  of  (l.l)  and  (1.2)  as  aany  tiaes  as  ve  need. 

In  classical  Mechanics  the  paraxeter  t  Is  identified  as  ties. 
Here  it  is  considered  as  the  fourth  coordinate  in  space -tine 
x*1  «j  (xk  ,  ct)  ,  (sst  c  =  l)  which  has  the  Metric 


d»  /v  = 


10  0  0 
0  10  0 

0  0  10 

0  0  0  -1 


(c  »  1) 


L.  J 

tbb  following  convention  is  used:  The  snail  Latin  subscripts  or 
superscripts  will  always  assuae  the  values  1,  2,  3  .  They  will 
signify  the  spatial  coordinates  of  the  space -ties  of  events.  The 
snail  Orttk  subscripts  or  superscripts  will  always  assuae  the  Values 
1,  2,  3,  k  .  They  will  designate  the  coordinates  of  space -tine. 

They  will  be  raised  and  lowered  by  the  aetrlc  .  The  large  Latin 

subscripts  or  superscripts  as  suae  the  values  1,  2,  3  and  will  be 
raised  and  lowered  by  the  Krone cker  delta  5^  .  Jhey  will  signify 
the  coordinates  of  the  reference  state. 

Slallikr  to  (1.2)  the  inverse  notion  of  a  continuous  aedlunaay 
be  described  by  three  functions 


(1.5) 


x’V)  (1.1,2,3)11.1,2,  5,  *) 


The  domain  of  X*()  Is  the  material  tube  swept  out  by  the  body: 

D  ■  (x^ :  c  B  ,  0  <  t  <  •)  in  space  -  time.  The  range  of 

3T()  Is  the  set  (B).  In  the  following  It  will  be  implicitly  assumed 
that  £  c  D  and  £  <  B  and  no  further  reference  will  be  made  to 
this  fact. 

The  three  functions  A)  are  assumed  to  be  invariant  functions 
of  x*1  under  the  group  of  homogeneous  Lorentz  transformations ,  £  , 


(1.6) 


*  -  aVv 


A01  A  V  -  A  a  AV 

A  v  %  Av  A  0 


,a 


that  is: 

(1.7)  Ax11)  -  **(**) 

Mathematically  this  is  all  that  is  needed  to  develop  a  rela¬ 
tivist  ically  invariant  kinematics.  However,  for  a  physical  theory, 
some  meaning  must  be  attached  to  the  functions  **().  Physically, 
consider  an  undeformed  body  B  which  before  deformation  is  at  rest 
in  some  frame.  In  ;hat  frame  the  body  B  is  described  by  a  set  (B) 
in  the  three  dimensional  space  with  coordinates  Xs  .  The  three 
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World  Li nea.  Velocity,  Deformation  Gradients 

It  is  now  meaningful  to  define  the  world  line  of  a  particle  A1 
by  the  curve  in  space  -  time  defined  by 

(1.9)  AK  •  X^x11) 

Ihe  particle  velocity  is  defined  as 


(1.10)  vk  s  fo— -r'"-' 


dx 


The  world  velocity  vector  is 
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functions  X^()  describe  the  mapping  of  the  set  D  onto  B  .  It 
will  be  assumed  that  such  a  reference  state  exists.  Thus  X*  can 
be  used  to  describe  lengths  and  angles  in  that  frame. 

It  is  also  assumed  that  (1.3)  i®  invertible  in  the  first  three 
variables  x  ,  that  is 

k  k/vK 

(1.8)  det  (xk^K)  /  0 

k  .  x^) 

X’K  '  ‘  ^ 


I 


1* 


Thus  the  following  useful  Identities  are  derived; 


(1.1*) 


S  & 


k 

i 


The  last  equation  yields  the  identity 

(1.15)  ua  x1^  s  o 


Material  Surfaces  and  Volume a 

It  is  convenient  to  introduce  the  invariant  derivative 


(1.16)  D  ♦  5  ua  ♦  a 


If  ♦  is  a  tensor  under  Lorentz  trana  format  ions ,  Dt  is  a  tensor 
under  Lorentz  transformations.  This  is  the  relativistic  generaliza¬ 
tion  of  the  material  derivative  ~  ,  in  fact 


The  concept  of  a  material  surface  is  important  for  the  formu¬ 


lation  of  the  dynamical  lavs  of  a  continuous  medium. 


■«*$ 
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t 

Def .  A  surface  f(xa)  =0  is  a  material  surface  if  uti  only  if 

(1.17)  f(«k(xK,  x4),  x4)  s  rex*) 

A  necessary  and  stiff icient  condition  for  the  surface  f (x°)  *  0 
to  he  material  Is  that 

(1.18)  Df  ■  0 

The  proof  is  obvious. 

Prom  (1.17)  or  (l.lB)  it  follows  that  if  the  surface  determined 
by  f(x°)  =0  is  a  material  surface,  the  surface  f(xa)  =  0  ,  where 
xa  is  related  to  x°  by  (1.6),  is  also  a  material  surface.  That  is, 
the  concept  of  a  material  surface  is  invariant  voder  Lorentz  trans¬ 
formations  . 

A  material  volume  is  a  three  dimensional  volume  containing 
material  points. 


Decomposition  of  Vectors 

A  four  vector  f°  is  called  space -like  if  f°  f Q  >  0  ,  null  if 
f'  ?  *  0  ,  and  time -like  if  f Q  f a  <  0  .  Prom  (1.11)  the  four  vector 
ua  is  a  unit  time -like  vector.  In  the  sequel  It  will  be  found  convenient 
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to  decompose  all  vectors  and  tenaora  Into  apace -like  and  time -like 
component a . 

An  arbitrary  four  vector  can  be  decomposed  into  a  time -like 
component  parallel  to  u3  and  a  apace -like  component  perpendicular 
to  ua  .  To  thia  end  ve  introduce  the  projectors  8° 


(1.19) 


Sap  *  6ap  t 


a 


err 


a 


,a 


V  5  V7  8t  =  6p  +  u  up  ■  “P 


The  projectors  aatiafy  the  identitiea 


(1.20) 


SQ  8P  =  8° 

P  7  ? 


SP_.  -  uP  8®  =  0 


In  genera],  the  four  vector  F°  can  be  uniquely  written  in  the  form 

I**  .  f“  +  uaf 

(1.21) 

f°u  »0 

a 


From  (1.20)  and  (l.2l)  it  follcva  that 


(1.22)  f  «  -  uQ 

(1.23)  f°  -  Sap  ^ 


It  can  be  eaaily  seen  that 
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***«  *  0 


That  is  f  is  space -like  or  a  null  vector  (light  signal). 
From  (1.15)  and  (1.19)  it  follows  that 


(1.210  --  * 


Therefore  the  vectors 


are  space -like, 


Deformation  Tensors 


The  deformation  of  bodies  can  he  described  by  the  deformation 
gradients  X*  .  The  invariant  strain  measure  (Toupin  [ll] ) 


-1 


KL 


(1.23) 


=  ,*** 


»<*  ,p 


c1kl  =  c“ 


is  the  relativistic  generalization  of  the  inverse  Grew  deformation 
tensor.  These  are  six  Lorentz  invariant  scalars.  From  (l.lV)^ 

(1.26)  ^  •(#«  -  v1  yi)  Xs  ,  X1  , 

-V 

From  (1.26)  it  is  clear  that  C  is  the  inverse  Green  deformation 
tensor  of  a  local  instantaneous  rest  frame. 


f 


l8 


From  (1.26)  vq  calculate 

Al 


(1.27) 


det  (C^)  =  det  (61J  -  v1  v^)  [det  (X*  .)Y 

=  (1  -  v2)  [det  (X*  .  )]2  >  0 


The  inequality  foil  ova  from  (l.'8  )g  and  the  physical  assumption  that 

g 

no  body  can  move  faster  than  the  speed  of  11 ght  (v  <  l)  .  The 


matrix 


Is  Invertible  and  It  is  easily  verified  that 


ckl  1  <6ij  ♦  «i  V  x\: 


-1  ITT 

is  the  inverse  of  C  ^  ,  i.e., 


(1.29)  CkmCML  =  >Cmc  *  BLk 


The  six  quantities  are  scalar  invariants  under  the  Lorents 

~  iCL 

group  since  they  are  invariant  functions  of  the  invariants  C 

The  matrix  has  a  very  simple  physical  meaning:  it  gives 

the  changes  in  length  and  angle  due  to  the  deformation  as  vleved  hy 
an  observer  in  a  local  instantaneous  rest  frame.  To  an  observer  in  a 
local  instantaneous  rest  frame  v*  »  0  ,  u*  *  0  and  (1.26)  and  (l.28) 
revert  to  the  forms  known  to  us  from  non-re la tirris tic  continuum 
mechanics  for  which  the  connection  of  to  length  and  angle 

changes  is  well  known  (cf.  Eringen  [  7  ,  Art.  7]). 


ika.nu  i*.JJUJJ!u,  Attl  WlMf 
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Theorem: 


The  quantities  (K  fixed) 


are  four  Teeters. 


Proof: 


AO 

x 


Prom  (1.34)  we  have 


a 

K 


r 


Since  Saa  is  a  tenaor 
p 


n  ■  ^  v  s76 


Thus 


a  6  7  JC  .  $a  a  7  “a 

A7  P  K^,6'3P  '  A  ~  'r  “ 


7  K  ,p  P 


since 


*.  ■  V 


But  from  (1.34), 


Consequently 


ia  y 


Aa  x'  XK  . 


P  L  7  K  ,P  L 


A°  xX 
7  L 
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or  in  view  of  (1.3*0^ 

K  A  P  *  K 

which  is  the  proof.  Thus  is  a  tvo  point  tensor.  (The  fact 

that  x^K  is  a  vector  under  change  of  X^,  with  p  fixed,  is  more 
or  less  obvious). 

Ike  relativistic  generalizations  of  the  Cauchy  def oration 


tensors 

are: 

(1.55) 

C<#  5 

V  *L,, 

(1.36) 

HI 

A 

»n  \  A  •  ° 

where 

iu  - 

cKL  k  l 

0  x  „  X  — 

,a  * 

From  these  and  (1.3*0  we  can  establish  the  identity 
(1.37)  F  *  S“p 

The  following  interpretation  can  be  given  to  *  An  infinites¬ 
imal  spatial  measurement  by  an  observer  in  a  Lorentz  frame  is  the 

v’  4  k  k  h 

separation  of  two  simultaneous  events,  (x  ,  x  )  and  (x  +  dx  ,  x  )  . 
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Since  tvo  events  simultaneous  in  one  Lorentz  frame  are  not  necessarily 
simultaneous  in  another  Lorentz  frame,  the  act  of  spatial  measurement 
is  not  an  invariant  operation.  The  difference  between  tvo  events 
x^  and  x^  +  dx^  is  a  spatial  msasuresMnt  in  the  local  instantaneous 
rest  fraae  at  x^  if  and  only  if 

Up  dx^  ■  0  at  £ 


If  the  tvo  pairs  of  events  [x^  ,  +  d?]  and  [x^,  x^  +  dx^]  are 
spatial  measurements  of  tvo  material  lines  in  the  local  instantaneous 
rest  frame  at  &  ,  then 


is  the  cosine  between  the  material  lines  before  deformation. 


Invariants  of  Strain,  Volume  Changes 

By  use  of  the  defining  equations  (1.33)*  (1*36)  of  the  Cauchy 
defor-jetion  tensors  and  (1.25)  and  (1.32)  of  Green  deformation  measures 
the  relations  of  invariants  cf  £"  and  £  to  those  of  material  strain 
measures  can  be  deduced.  Thus, 


'<r, 
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(1.38) 

tr  £  * 

tr  £  , 

tr  £2 

:  tr<*)2 

,  tr  £5 

=  trfA5 

(1.59) 

tr  £  • 

-1 

tr  £  , 

tr  / 

«  tr  (C)2 

>  tr  £5 

■J-  * 

■  tr  (C)* 

where  we  used  the 

notation 

tr  t 

■*. 

,  tr  £ 

«  ca 
a 

2 

,  tr  £ 

m 

> 

«> 

r  •  •  • 

-1 
tr  £ 

* 

K 

,  tr  £ 

*  cK 

C  K 

»  •  •  • 

Suppoee  that  the  material  volume  dV  is  deformed  into  dv  in 
a 

the  I/srentz  frame  x  .  In  an  instantaneous  rest  frame  coinciding 
vith  xa  the  deformed  volume  is  dvQ  given  by 


(1^0) 


dvo 


where  J  is  the  Jacobian  defined  "by 


(l.M)  J  -  £«oe76  «kx«  ^.p  *,r  "6 


Alternatively 


j— 

0 
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where 


(l-W  *  -  5  «*“  v  X°K  *7M  "6  “  I 


From  (1.4l)  and  (1.43)  ye  can  also  deduce  the  Identities 


(1.44)  u6  «  J  (vfu  Xs  _  X1*  „  X14 


O07& 


KLM  “  ,a  “  ,7 


t1-^)  ‘a*  “  J  •<**  x“k  A,  A  * 


The  above  Identities  are  useful  In  the  reduction  of  constitutive 
equations  for  Isotropic  materials. 


Compatibility  Conditions 

In  classical  continuum  mechanics  the  compatibility  conditions 
for  the  deformations  tensors  are  deduced  from  the  requirement  that 
the  space  remain  Euclidean.  This  is  not  the  case  in  the  relativistic 
theory  presented  here.  Suppose  the  deformation  tensor  c^  is  given 
at  each  event  (&)  ,  then  the  system  of  equations  for  X^(jt) 

°ce  -  6klxVL,p 

is  over  determined.  Differentiate  (1.46)  with  respect  to  x7  . 

' 
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(1>7)  eop,r  -  SKL  *„  *,»  +  6kl  **  «  ^.Pr 


Define  a  Chris toff el  symbol  of  the  second  kind  by 


(1.48)  IflP.y]  *  |  +  %,a  -  'op.r1 


Then  using  (1.47),  (1.48)  ve  see  that 


(1.49)  [oe,7]s  -  hL*,ce*,r 

But  (1.49)  is  equivalent  to 


(1.50)  X*  -  Xs  (  7  ) 

,ap  ,y  up 

where 

(1.51)  (a7p)  i  s78  [ce,8)£ 

Is  called  Christoffel's  symbol  of  the  second  kind. 

The  system  (I.50)  is  completely  integrable  if  and  only  if 


(1.52) 


*»,»  *  f.7  '.V  <=%> 

■  *>,.  *  1.V  <«%> 


(See  Eisenhart  [29],  p.  1  and  pp.  186-188. ) 
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Thu is  a  necessary  and  tuff  1  clast  condition  that  (1.50)  la  completely 
integrable  la  that  the  "curvature"  tana  or  vanish,  i.e. 


(1.53)  ^  *  0 

vhara 


<!•*>  W  8at[(PV,r-(PV,8+(<rV(PV 


'aV'P°7)] 


This  la  the  relativistic  generalization  of  the  compatibility  condi¬ 
tion*  for  . 

cp 


Rates  of  Deformation  Gradients  and  Strain 
The  following  lens*  la  useful  in  the  kinematics  of  continuous 

media 


(1.53) 


To  prove  this  we  note  from  the  definition  (l.l6)  of  the  operator  D  that 


(»“  Xs  ) 

ta  tp 


I JC 


According  to  (1.15)  the  first  term  on  the  right  hand  side  is  zero  and 
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we  get  (1.55 )• 

A  dual  to  (1.55)  i« 

(1.56)  Sap  D  xPK  ■  xeK  u°jP 
For  the  proof  of  this  we  use  (1.34  )g. 

d^/l1  ■  D  x\  +  r\e  L  x\  -  0 

Hiltiplying  this  by  xaR  and  using  (1.55)  end  (1.34)  we  obtain  (1.56). 
It  is  useful  to  introduce 

U.57)  ^  -  Sg 

in  terns  of  which  (1.56)  n»y  be  expressed  as 
(I.5«)  s“p  D  x\  .  S“P  xeK 

Using  (1.34)^  we  have 

n  a  a  _ 

ua  D  x  K  -  -  x  K  D  ua 

When  this  and  the  expression  (I.19)1are  employed  in  (1.58)  we  find 

■  ^  4  +  D  U@ 


(1.59)  D  x“K 
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The  invariant  derivative  of  nay  now  he  calculated  by 

(1.60)  DCia  -  a^*“K^L 

where 

(1-6D  %  ■  I  <  V  .  5^, 

The  proof  of  (1.60)  is  immediate  from  taking  the  invariant  deriva¬ 
tive  of  (1.32)  and  using  (l.59)» 

* 

The  tensor  4  is  the  relativistic  generalization  of  the  deforma¬ 
tion  rate  tensor,  It  is  clear  that  for  locally  rigid  motions  D  -  0  ; 

consequently,  we  have 

Theorem:  A  necessary  and  sufficient  condition  for  locally  rigid 

motion  is 

(1.62)  5^  -  0 

This  is  the  relativistic  generalization  of  the  Killing's  theorem  of 
differential  geomstry. 

The  relativistic  generalization  of  spin  is  defined  by 
d.6j)  2^  «  \  (Sqp  -  u^)  -  S[Qf5) 


Upon  adding  thia  to  (l.6l)  va  have 


d.(»)  +  2^ 

The  invariant  derivative a  of  various  a train  measures  and  other 
tensors  can  be  calculated  by  use  of  the  apparatus  aet  above.  Here 
va  tflve 


(1.63) 


®  CQp  +  “ 


,P  ear 


+  u' 


,a  erP 


0 


which  la  obtained  by  taking  the  invariant  derivative  of  (1.??)  and 
using  (1.55)-  Eii»  expression  will  be  found  useful  in  the  treatment 
of  isotropic  materials. 

Another  useful  result  is  the  identity 

(1.66)  Dj  -  -  J  uaa 

which  can  be  proven  by  (1.43)  and  (1.59). 
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RELATIVISTIC  BALANCE  LAWS  FOR  A  COMTOJOCfS  MEDIUM 


General  Balance  Lavs 


i 

In  nonrelativistic  continuum  mechanics  the  balance  lavs  are 
written  in  the  form 


(2.1) 


ST  /  *dv 

T(t) 


■  /  £•<!£+/  s  dv 

*(t)  v(t) 


where  $  is  some  quantity  whose  influx  is  £  and  whose  supply  is  s  . 
Equation  (2.1)  holds  for  an  arbitrary  material  volume  v(t)  whose 
enclosing  surface  is  ^(t)  .  This  equation  can  be  integrated  from 
ti  to  tg  for  arbitrary  t^  <  tg  to  obtain 


(2.2) 


Under  the  proper  smoothness  conditions  the  left  hand  side  of  (2.2) 
can  be  replaced  by  an  integral  over  the  three  dimensional  circuit  in 
space -time  enclosing  the  four  dimensional  volume  (called  material 
tube)  swept  out  by  the  material  volume  v(t)  in  the  interval 

(*1'  *2^  •  Hcnce 


/  *“  *•*  *  /'  dr* 
dv^  ■  dv  dx1* 


(2.J) 
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where 

♦*  .  -  ik  +  *  vk  (k  •  1,  2,  5) 

(2A) 


and  d#^  and  dv^  respectively  denote  the  elements  of  three  diuan- 
sional  oriented  surface  and  the  volume. 

To  obtain  (2.3)  we  recall  that  a  material  volume  in  space -time 
is  a  three  dimensional  surface  with  peramet;  izatlon 


(2.5) 


K  ■  1,  2,  3 

(t  =  constant) 


Also  a  three  dimensional  surface  in  space -time  swept  out  by  a  material 
surface  can  be  parametrized  by 


(2.6)  xk  »  xk  (Xs  (ux,  u2),  u^)  ,  =  u5 


so  that  the  oriented  three-dimensional  surface  element 
expressed  by 


can  be 


(2.7) 


ds 


3u 


a  c 


[a.  0 

dx  3x 
V&7  Su^  ^ 


7  1 

^  du1  du2  du^ 


where  a  bracket  enclosing  indices  indicates  skew-symmetry  as  defined  by 
a[0P7J  .  1  (A^  +  +  APra  .  aP07  .  A7?a  _  a07P) 
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m 

The  quantity  ♦  will  be  said  to  be  conserved  if  i*  :  0  ,  s  =  0  . 
From  (2.3)  and  (2.4)  the  quantity  ♦  is  conserved  if  and  only  if 

/*“**  3«  '  ° 

(2.8)  -  *ou* 

♦  «  -  «  ♦  Vl  -  v^ 

o  a 

So  far  all  this  holds  for  either  classical  continuum  mechanics  or 
relativistic  continuum  mechanics.  lhe  only  distinction  betveen  them 
is  the  group  of  transformations  in  space -time  for  which  the  integral 

(2.9)  •  [♦,  »51  -  J ♦ads^a 

C3 

has  a  specified  transformation  law  —  t^e  Galilean  group  for  classical 
mechanics  or  the  Lorentz  group  for  relativistic  mechanics.  In  the 
following  subsections  the  lavs  of  conservation  of  mass,  balance  of 
ener.'^y -momentum,  balance  of  angular  momentum,  and  the  second  law 
of  thermodynamics  are  formulated  in  a  relativistically  invariant 
manner.  In  general,  these  laws  will  have  the  form  (2.3)  with  a 
transformation  law  for  a  set  of  quantities  of  the  form  (2.9)  specified 
for  the  Lorentz  group. 


\ 
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'  Conservation  of  Particle  Kuaiber 

In  nonrelativistic  mechanics  the  postulate  of  conservation 
of  mass  and  the  conservation  of  the  number  of  particles  are  equivalent 
for  nonreacting  systems.  It  is  veil  known  that  this  is  not  th*  case 
in  relativistic  mechanics.  In  relativity  the  mass  is  closely  related 
to  the  energy  of  a  particle.  In  general  the  mass,  even  the  rest 
mass,  of  a  particle  varies.  The  concept  of  the  non-creation  or 
indestructability  of  particles  still  remains  valid  in  classical 
relativistic  theories  of  nonreacting  substances  (the  word  classical  is 
used  in  opposition  to  relativistic  quantum  theories  where  the  idea 
of  creating  or  destroying  particles  is  essential).  It  is  assumed  that 
the  number  of  particles  contained  in  a  material  volume  V(t)  is 
a  constant  of  the  motion. 

To  formulate*  the  law  of  conservation  of  particle  number  we 
assign,  to  every  three  dimensional  surfac  :  s^  ,  a  positive  scalar 
H[s^3  of  the  form 


* 

Each  lav  formulated  in  these  subsections  is  deduced  along  the 
lines  leading  from  (2.1)  to  (2.3)  from  the  corresponding  classical 
law.  Though  the  correspondence  (2.4)  is  usually  listed,  it  should  be 
noted  that  this  is  not  needed  for  a  relativistic  formulation  of  the 
laws  of  mechanics  ,  Only  the  form  of  the  equations  is  necessary.  For 
a  physical  feeling  and  interpretation  of  the  various  components  of 
space -time  tensors  and  vectors  the  correspondence  (2.4)  is  indispensable 
and  will  usually  be  given. 


5^ 


(2.10)  H[«3)  =  J  na 

*3 

which  is  postulated  to  be  Invariant  under  the  group  of  Lo rente 
trann format ions.  The  lav  of  conservation  of  particle -number  states 
that  for  every  material  volume,  V(t)  ,  Nls^]  is  conserved. 

From  (2.8)  this  is  equivalent  to 

J  na  ds^  ■  0 

(2.11 )  na  =  nQ  ua 

n  ■  -  na  u  >  0 
o  a 

where  the  integral  is  over  an  arbitrary  material  tube.  The  scalar 
nQ  is  called  the  rest  frame  particle  number  and  is  related  to  the 
particle  number  ,  n  ,  by 


(2.12)  n 


The  use  of  Green -Gauss  theorem  in  the  integral  equation  (2.11).  leads 
in  the  usual  way,  to  the  following  differential  equation  and  Jump 
condition: 


(2.15)  n°  „  «  0  ,  [na]  across  E  (x^)  *  0 

where  [f]  ■  f  +  f”  denotes  the  Ju mp  of  f  at  any  discontinuity 
surface  Z(x^)  »  0  . 


55 


It  is  convenient  to  write  (2.1J)  u 

(2.14)  0»o*»oJBp  =  0 

where  the  following  identity  ha*  been  employed 

(2.15)  SPp  *'  uPp 

By  using  (1.66)  one  can  show  that  (2.1*0  has  a  solution  of  the  form 

(2.16)  n0  =  n°J 

where  n°  is  a  function  of  X*  . 
o 

In  the  following,  various  scalars  ♦  will  appear  in  the  balance 
equations.  It  is  convenient  to  define  another  scalar  #o  related 
to  ♦  by 

(2.17)  ♦  --  no«o 
From  (2.1*0 

(2.18)  (♦up)p  =  n0W0 

This  identity  is  used  frequently  in  the  formulation  of  the  remaining 


balance  lavs. 
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Balance  of  Energy- Momentum 

In  relativity  theory  the  lav?  of  balance  of  momentum  and  energy 
are  closely  connected.  They  are  the  components  of  a  four  vector. 

For  physical  reasons  It  Is  expected  that  the  time  rate  of  change  cf 
the  momentum  Is  equal  to  the  forces  applied  to  the  body  and  that 
the  time  rate  of  change  of  the  energy  Is  equal  to  the  vork  done  on 
the  tody.  In  a  classical  nonrelatlvlstlc  theory  of  a  continuous 
medium  these  lavs  take  the  form:  the  time  rate  of  change  of  the  mo¬ 
mentum  contained  In  a  material  volume  Is  equal  to  the  sum  of  the 
forces  applied  to  the  volume.  These  forces  can  be  decomposed  Into 
tvo  parts,  one  arising  from  surface  tractions  and  the  other  from  the 
body  forces.  The  time  rate  of  change  of  the  energy  contained  In  a 
material  volume  Is  squal  to  the  heat  flov  through  the  surface  plus 
a  heat  supply  Inside  the  volume  plus  the  vork  done  by  the  surface 
tractions  and  body  forces.  In  nonrelatlvlstlc  mechanics  the  momentum 
Is  usually  due  to  the  motlcn  (kinetic  momentum)  and  the  snergy  Is 
the  sum  of  the  kinetic  energy  and  the  Internal  energy.  In  relativistic 
theories  It  is  possible  to  have  momentum  of  nonkinetic  origin.  There 
Is  no  advantage  In  decomposing,  a  priori,  the  snergy  momentum  tensor 
Into  kinetic,  thermodynamic,  and  other  parts. 

The  energy -momentum  of  a  material  body  is  determined  by  assigning 
to  every  three-dimensional  subspace,  s^  ,  of  a  four -dimensional 
material  tube  four  functions  F^[s^3  vhich  are  components  of  a  four 
vector  under  Lorentz  transf ormations .  The  functions  have 
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the  fora: 


(2.19)  *  J  l*‘v  <ujv 

*5 

vhere  T*iV  is  a  second  order  tenser.  The  balance  of  energy -momentum 
states  that  for  every  material  tube  we  have 

(2.20)  j  ^  d»jv  *  J  1?  dv4 


where  the  four  vector  f^  represents  the  body  forces  and  energy 
supply  per  unit  volume. 

In  nonrelatlvlftlc  continuum  mechanics  we  have  the  Identifications: 


(2.21)  T14  =  p1  ,  T1^  ■  -t1J  +  p1  vJ  ,  T44  «  e  , 

Ju  .  1  +J1  .  .  1 

i  -  q  -  v  Vj  t  e  v 


and 


(2.22)  f1  «  f1  ,  t 4  =  h  +  t  •  z 

vhere  p1  Is  the  momentum  density,  t1^  Is  the  stress  tensor,  e  Is 
the  energy  density  (internal  plus  kinetic)  and  q1  is  the  heat  flux, 
f1  is  the  body  force  and  h  is  the  body  heat  supply.  Ibis  identifi- 
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cation  it  useful  for  establishing  a  familiarity  with  the  physical 
meaning  of  the  components  of  the  energy-momentum  tensor j  however, 
it  is  not  essential  for  a  treatment  of  relativistic  mechanics.  The 
only  assumption  employed  in  decomposing  the  energy -skomentum  tensor 
is  the  existence  of  a  world  velocity  vector  in  the  four  dimensional 
material  tube  in  space -time. 

In  the  usual  manner  (2.20)  leads  to 


(2.25) 


*  0  across  Ztx*1)  :  0 

•w  s  V 


The  tensor 


can  be  decomposed  into  a  scalar,  two  spatial 


vectors  and  a  spatial  tensor  by  applying  the  projection  operators 
g<* 

3  fi  * 


(2.24) 


this 

end  define 

w 

=  ^uau. 

a 

<1 

t78 

a 

P 

& 

II 

1 

3^ 

0 

u 


It  is  easily  shown  that 

(2.25)  1*  =  o.uauP  +  ua,P  +  pauP  -t* 
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It  should  be  noted  that  this  decomposition  is  perfectly  general  and 

depends  only  on  the  character  of  the  world  velocity  field.  We  have 

a  B 

the  following  physical  interpretations :  The  tensor  w  u  u  is  the 

kinetic  energy -momentum  tensor  with  the  mass  density  given  by  the 

2 

famous  formula  of  Einstein  E  =  me  .  Thus  there  is  a  contribution 
to  the  mass  density  due  to  the  internal  energy  of  the  body.  The 

Q 

four  vector  qp  in  the  local  Instantaneous  rest  frame  reduces  to 

i  i  B 

[q  ,  0]  ,  where  q  in  the  heat  flow  vector.  Thus  qp  is  called 

the  heat  flow  four  vector.  The  four  vector  pa  becomes  (p1,  0] 

in  the  local  Instantaneous  rest  frame,  where  p*  is  the  nonmechanical 

momentum.  Thus  p°  is  called  the  nonmechanical  momentum  four  vector. 
06 

The  four  tensor  t  is  the  relativistic  stress  tensor  since  it 
reduces  to 

o' 

0  0 

in  the  local  Instantaneous  rest  frame. 

From  (2,24)  and  (l.20)  the  following  identities  follow: 


(2.26) 


a 
<1  u« 


a 


a 

P  u, 


a 


0 

0 


0 


The  fourth  component  of  (2.23)^  is 
quently  the  first  law  of  thermodynamics 


the  balance  of  energy.  Conse- 
i8  contained  in  these  equations; 


HO 


however,  it  i8  not  the  fourth  component  of  (2.23).  To  obtain  the 
first  law  of  thermodynamics  one  nuit  take  the  projection  of  (2.23) 
onto  ua  .  This  operation  in  the  relativistic  generalisation  is 
the  counterpart  of  taking  the  scalar  product  of  Cauchy's  equations 
and  the  velocity  and  subtracting  the  result  from  the  energy 
equation.  We  recall  that  this  operation  eliminates  the  kinetic 
energy.  Hence 


*****  *  *** 


or 


(2.27)  (ua  1»>  -  T*3  ua<(5  -  f“  ua 


But  from  (2.25),  using  (2.26),  we  have 


(2.28) 


ul®  .  .  .  ue  - 


*  •  »*  %  •  *® 


Ihus  equation  (2.27)  reduces  to 


(2.29) 


-(u  uP) 


-  q 


a 

-  p  Du 


a  +  ^  ua,p 


-  f 


a 


Define  c  by 


■  no* 


(2.50)  U 


From  (2.l4),  (2.29)  becomes 


(2.31) 


no  ^  +  Ae +  p£ 


%  *  ^  ua,p  "  *  f\ 


IMs  is  the  relativistic  generalization  of  the  first  lav  of  thermo¬ 
dynamics  for  a  continuous  medium  (see  Eckart  [ 17 ] ) . 

The  generalization  of  Cauchy's  lave  can  be  obtained  by  applying 
S7a  to  (2.23).  Ihe  result  is: 


Only  three  of  the  four  equations  (2.32)  are  independent. 


Principle  of  Moment  of  Energy -Mome ntum 

In  the  classical  theory  of  continuum  mechanics,  for  nonpolar 
materials,  there  is  a  balance  lav  of  angular  momentum  vhich  by  the 
arguments  presented  at  the  beginning  of  this  chapter  can  be  written 
in  the  form? 

(2.33)  j  x[j  T110^  *  J  f11  dv4  ,  (i,J  =  1,  2,  3) 

for  an  arbitrary  material  tube.  (2.33)  is  the  spatial  part  of  the 
four  tensor 


(2. 5*0 


li-2 


j  *[V1(1  da^  =JVa  dv4 

If  it  Is  required  that  the  angular  momentum  is  balanced  in  every 
Lorentz  frame  then  (2.33)  implies  (2.3*0.  She  remaining  three  com¬ 
ponents  of  (2.3*0,  other  than  (2.33),  i«e., 

(2.35)  j  *[k  T1^  ds^  «  J  ^  flJ  dv^ 

are  the  expressions  of  equivalence  of  energy  flux  and  momentum  flow. 

If  (2.35)  is  assumed  to  hold  in  every  Lorentz  frame,  (2.33)  and  (2.3*0 
must  hold.  This  would  indicate  that  (2.3*0  must  be  a  basic  law  for 
nonpolar  materials. 

That  there  should  be  a  shew -symmetric  four  tensor  balance  law 
to  replace  (2.33)  is  indicated  also  by  the  following  argument: 

Consider  a  closed  conservative  system.  In  modern  physics  the  con¬ 
servation  of  momentum  is  interpreted  as  the  invariance  under  spatial 
translation,  while  the  conservation  of  energy  is  implied  by  the 
invariance  under  time  displacements.  In  classical  mechanics  this 
leads  to  one  vector  law  and  one  scalar  law  because  the  Galilean  group 
is  used.  In  relativistic  mechanics  this  leads  to  a  four  vector  equa¬ 
tion  since  the  Lorentz  group  is  used.  In  classical  mechanics  the 
Galilean  group  allows  spatial  rotations  and  invariance  under  these 
leads  to  the  conservation  of  angular  momentum.  The  Lorentz  group, 
however,  allows  rotations  in  the  four  dimensional  space  of  events  and 


thus  Implies  three  extra  conservation  lavs.  (See  Bogoliubov  and 
Shirkov  [21]  for  an  excellent  discussion  of  these  points.)  Con¬ 
tinuum  mechanics  treats  open  systems.  The  conservation  lavs  of  a 
closed  system  are  replaced  by  corresponding  balance  lavs.  The 
balance  of  energy -momentum  for  a  continuous  medium  has  already  been 
formulated.  The  following  assumption  formulates  the  principle  of 
moment  of  energy -momentum  which  serves  as  the  relativistic  generali¬ 
zation  of  the  balance  of  angular  momentum. 


To  every  three  dimensional  subspace,  s^  ,  of  the  four  dimensional 
material  tube  assign  a  skew  symmetric  tensor  function  if*  1. 1 


(2.36)  b3 

M(°e)u  ,  o 


8°^  is  called  the  spin  tensor1  (cf .Papapetrou  [22]).  For  nonpolar 
Qfiu 

materials  s  ^  z  0  . 

1This  tecs or  includes  such  effects  as  the  couple  stress  and  the 
intrinsic  spin  of  continuum  mechanics. 


t  JiUWWjrSi»*  T~ - ir - 


•  WIH!iWBBBKWWtew8P«iM — uni—  m 


The  lav  of  balance  of  scssnt  of  energy -Decent  um  le  that  for 
every  material  tube 


f  if*  d.Si  =  JV°  f01  ♦  L<*)  dvt 


(2.38) 


L(Qe)  =  0 


where  l*3  is  the  four  dimensional  analogue  of  the  body  torque, 
Equation  {2.38)  lead*  to 


(a  jfiJ  +  Loe 


(2.39) 


=  x 

*(» 


[lf*]Z  «  0  acrosc  £(&)  *  0 


Using  (2.23)  and  (2.37),  (2.39)  l«ad*  to 


(2.U0)  3°^  -T[0P]  =  L*5 

For  nonpolar  materials  L01^  =  0  and  (2.U0)  reduces  to 


(2.4l)  T[ael  =  0 


From  (2.25),  (2.fcl )  is  equivalent  to 


(2.42)  t[QPl  .0  ,  p“  =  q“ 


4—  nil 
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Throughout  this  work  s^K  =  0  .  In  the  case  of  electromagnetic 

06 

interaction  with  a  material  "body,  a  body  torque  L  of  electro¬ 
ns. 

magnetic  origin  will  be  introduced.  Hie  inclusion  of  s^^  is 
necessary  if  one  wishes  to  formulate  a  relativistic  generalization 
of  mechanical  theories  of  couple  stresses.  A  fuller  investigation 
of  the  properties  of  this  tensor  is  left  for  further  research. 


Second  Lav  of  Thermodynamics 

In  modern  continuum  mechanics  the  second  lav  of  thermodynamics 
is  considered  to  be  a  restriction  on  the  form  of  the  constitutive 
equations.  However,  formulation  of  this  lav  is  independent  of  the 
character  of  the  media  under  consideration.  This  lav  is  expressed 
in  the  form  of  an  inequality  called  the  Clausius -Duhem  inequality 
which  has  the  form  (2.1)  with  the  equality  replaced  by  an  inequality. 

To  formulate  the  relativistic  extension  of  this  inequality, 
assign  to  every  three  dimensional  subspace,  s^  ,  of  a  material  tube 
a  scalar  invariant  H[s^]  of  the  form: 

(2.43)  H(s j]  «  J  »iadsj0 

S3 


The  second  lav  of  thermodynamics  is  the  statement 
(2.44)  J  1°di *5a+  f  r  dv^  £  0 


for  every  materiel  tube.  Here  the  scalar  Invariant  r  Is  the  supply 
of  entropy  from  extraneous  sources.  Inequality  (2.44)  leads  to  the 
local  lav 


(2A3) 


n  +  r  >  0 

,a 


[naJ  l  £  0  across  Z(xM)  s  0 

*•*  *•  .a 


It  is  convenient  to  decompose  Ti  into  its  space -like  and.  time -like 
component a . 


(2.46)  =  Sap  if  , 


n  =  -n  u„ 
o  a 


so  that 


(2.47)  »1M  =  tl0  uM  +  f 


A  simple  thermodynamics  process  is  one  for  which 


(2.48) 


M  <£  .  hO 

'  •  r  ’  r  =  - 


b  =  f  u  ,  e  >  0 
o  u 


The  quantity  6  is  a  scalar  invariant  called  the  temperature,  and 


“>*«■«««  j4*) 


h  is  the  heat  supply  term  Define  n  by 
o  oo 


(2.49)  n0  •  B0  r,^ 


*7 


lbs  second  law  (2.^5)  is  then 

(2.50)  n  Dn„  +  »P  _  +  h  i  0 

#00  >P 

For  a  simple  thermodynanic  process  (2.50)  he  cones : 

(2.51)  n0  Dn^  ♦  (£)>p  ♦  Ja.  *  0 

Ibis  is  the  form  of  the  second  lav  introduced  hy  Eckh&rt  [17]  for 
fluids.  Eliminating  hQ  from  (2. 31)  and  (2.51)  ve  get 

6  p°W  t°^  U 

(2.52)  v^oo-?^  -y,*--ir  +  -tm 


Die  inequality  (2.32)  is  useful  for  the  reduction  of  constitutive 
equations . 


CHAPTER  III 
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CONSTITUTIVE  THEORY  OF  MECHANICAL  MATERIALS 


General  Principles 

In  general  the  system  of  balance  lavs  proposed  In  Chapter  II 
is  inadequate  for  the  treatment  of  problems  of  relativistic  continue 
except  in  some  special  cases.  The  properties  of  the  medium  are 
brought  into  consideration  through  a  set  of  constitutive  equations. 

In  general  a  constitutive  theory  should  satisfy  certain  principles^: 

i)  Principle  of  Causality:  The  behavior  of  the  material  at 
the  event  x  is  determined  only  by  events  lying  in  the  past  light 
cone  at  £  •  That  is  only  by  those  events  x  which  satisfy  the 
following  inequalities. 

(x  -  i)  •  (x  -  &)  <0 

(&-£)•  Ji  £  0 

where  ip  the  world  velocity  at  x  . 

ii)  Principle  of  Locality:  The  behavior  of  the  material  at 
an  event  depends  strongly  on  the  properties  of  the  mater ialj^.  the 

neighborhood  of  the  event.  _  _ — *■"■ 

ill )  Lorentz  Invariance :  The  constitutive  equations  are  covariant 
under  the  orthochronous  proper  inhomogeneous  Lorentz  group,  i.e.,  the 


■Tor  a  discussion  on  non -relativistic  constitutive  theory  see 
Eringen  [7,  Ch.  V]  and  [9]. 
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group  for  which  det  £  =  +1  ,  A  ^  >  0  . 

iv)  Material  Invariance :  The  constitutive  equations  are  invariant 
under  the  symmetry  group  which  characterizes  the  material  in  the 
Legrangian  frame  . 

v)  Consistency:  The  constitutive  equations  must  he  consistent 
with  the  balance  laws  of  particle -number,  energy -momentum,  moment  of 
energy -momentum,  the  lave  of  electromagnetism  (to  he  formulated  in 
the  next  chapter),  and  the  eecond  lav  of  thermodynamics. 

vi)  Equipresence :  An  independent  variable  that  appears  in  one 
constitutive  equation  should  appear  in  all  constitutive  equations 
unless  excluded  by  one  of  the  above  principles  i.  -  v. 

A  word  ie  in  order  about  the  requirement  of  invariance  under 
Lorentz  transformations.  In  the  modern  theories  of  continuum 
mechanic e  the  constitutive  equations  are  covariant  under  rigid  body 
motions.  Although  it  appears  that  claes ically  the  invariant  group  of 
mechanics  is  the  Galilean  group  [50],  [51],  the  formulation  of  modern 
continuum  theories  decomposes  the  forces  acting  on  the  body  into 
external  forcee  and  internal  forcee  and  assumes  that  the  internal 
forces  are  objective  under  ri^id  motions  while  the  external  forcee 
are  not  objective  [52].  In  this  article  we  make  no  such  distinction 
between  forces.  The  objectivity  under  Lorentz  transformations  is 
adequate  for  the  theories  presented  in  the  following  sections. 

In  the  next  two  subsections,  relativistically  invariant  consti¬ 
tutive  theories  for  tbermoelaetic  solids  and  thermoviecous  fluids 
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are  presented.  These  are  staple  relativistic  generalizations  of  the 
corresponding  classical  theories  of  nonpolar  Materials  for  which 
«  0  ,  L°^  *  0  .  Thus  (2.42)  is  valid. 

Tfaeraoelastic  Solid 

For  the  construction  of  constitutive  equations  of  a  theraoelastic 
solid,  an  appropriate  set  of  independent  variables  is: 

(3.1)  e  ■,  Xs  ;  ep  ,  xp 

where 

(3.2)  ep  «  sp“  (e  a  +  e  &>a) 

* 

choice  of  6^  as  an  independent  variable  appears  to  be  a  natural 
one  through  the  examination  of  the  second  law  of  thermodynanics  (2.52) 
for  nonpolar  aaterlals: 

(3.3)  »0  <»%,-!»•>  -4  V*F  %  *  0 

V 

where  (2.42)  has  been  used. 

Ve  now  write  constitutive  equations  for  the  dependent  variables 

a  as 

«  ,  ,  q  ,  t  in  the  general  foras 
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where 

(5.8)  *Q  * 

(5-9)  Ox  »  xaK  ^ 

The  quantity  tQ  i»  called  the  free  energy. 

At  any  event  £  the  following  quantities 

(5.10)  e  ;  De  ;  X*  j  fl*  ;  D6K  ;  DX*^  (I  -  1,2,3  i  O  -  1,2, 3» 

are  Independent  in  the  sente  that  there  always  exists  a  action  such 

that  for  an  arbitrary  set  of  numbers 

(3.11)  T  }  B  J  l^a  ;  ^  ;  AK  ;  DKa  (a  -  1,2,3>  ;  l  -  1,2,3) 

subject  to  T  >  0  ,  det  (7°^  I*  FLJ  >  0  ,  there  exists  an  allowable 

a  p 

motion  such  that 

(3.12)  e  -  T  ,•  DS  -  B  ;  X*  a  -  /a  J  eK  -  T*  ;  D6K  -  AK  ;  -  DKa 

at  an  arbitrary  event  &  .  The  inequality  (3*7)  contains  terms  of 


the  form: 


(3-13)  g(  )  y  +  f(  )  >0 

vfaerc  y  ie  from  the  set  (3*10)  and  g(  )  and  f(  )  do  not  depend 
on  y  .  Since  y  is  independent  in  the  above  sense,  a  necessary 
and  sufficient  condition  that  (3.13)  holds,  for  arbitrary  values 
of  y  ,  is  that 

(3-lM  g(  )  «  0  ,  f{  )  0 


results 
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(5.i7)  t  =  nQ  -4-  S7(/ ,p) 

*** 

(5-l8)  dr  s7ic?,e]  z  0 

>  7 


The  most  general  form  of  the  constitutive  equations  for  an 


anisotropic  thermoelastic  solid 


Lng  the  requirements  of  Lorent? 


Invariance  and  nonnegative  entropy  production  is  the  following: 

i.  The  free  energy  lrQ  has  the  functional  form: 


(5.19)  ,0  ■  i“)  .  *■  ■  l"  „  I1 . 


ii .  The  entropy  and  the  stress  tensor  t**  are  determined 

from  the  free  energy  i|tq  by 

(5.20)  n  =  o 


K  T 

(3.21)  t  =  -2n  r  f  -2- 
op  0  ,a  ,p 

oC 


iii.  The  heat  flow  vector  q  has  the  form: 


(3.22)  ^ 

where 

(3.2^)  ^  *  Q  x(e  ,  clKL  ,  eK) 


~  i  .  i  |— ii|-TMinmrTi - mMiM 


I 
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and  satisfies  the  inequality 


(3.24)  ^  6*  £  0 


The  inequality  (5-24)  implies  that 


"2kl 


(3.25)  <^(6  ,  C1^  ,  0)  =  0 


The  proof  of  statement  (i)  follows  from  (3.13 L  and  Lorentz 


invar  ic.nce.  According  to  (3.15 )g  have 


*0  *  *o(e  ’  **.(,> 


Nov  tQ  is  an  invariant  function  so  that 


V9  *  A«  *V  •  V9'^ 


.<x 


for  an  arbitrary  Lorentz  transformation  A  .  There  exist* 

P 


at  each  event  jc  a  Lorentz  transformation,  &  ,  such  that 


i 


5  X* 


4 


=  0 


,0 


From  (l.lfc)  the  matrix  .  is  invertible,  therefore  it  has  a 
unique  polar  dec  exposition 


.i  9  ^ 


\ 

> 


where 


yKL  yM  ,  -JKM 

jj 


and 


"u 


"u 


R*1  =  fK 


'W 


^  --  *i 


Thus 


RUA/xKfl  .  v® 

1  JP 


Therefore  at  each  event  jc  ,  there  exists  a  Lorentz  transformation 


A  -  such  that 
P 


?®  .  h^A/  ,  A1*  .  ^ 


and 


a-x*,  *  v1*  ,  3*i«  *  o 


,P 


,P 


Thus  t  can  be  considered  as  a  function 
o 


♦o  1  *0(e  *  vKL> 

JCL  -ivT 

Since  v  is  a  function  of  C  we  have  the  proof  of  (i).  The 
proof  of  the  first  part  of  (ii),  (eq.  3.20),  was  already  given  by 
(3.1$).  Hie  second  part  (eq.  3.21)  follows  using 


* 


—  (Xs  Q  &Lm  +  X1 


,P  M  M 


-  .  - 


in  equation  (5. IT).  Equation*  (3»l6)  and  (3.l8)  are  now  satisfied 
identically.  The  proof  of  statement  (iii)  follows  from  the  fact 
that  0^(6  ,  p  ,  eK)  is  an  invariant  function  under  Loren tz 
transformations : 


V6  '  Aa  xK,p  '  ®K)  =  Ve  '  xK,a  '  6*> 


Using  the  same  argument  that  was  employed  above  for  must 

have  the  form: 

Qj  =  ^(e,  c1kl  ,  eK) 

Therefore  (5.22)  and  (5*25)  follow. 

The  inequality  (5.2^)  implies  that  if  ©2  =  =  0  , 

>  0  if  <  0  and  Q1  <  0  if  6^  >  0  .  If  we  assume  that 
JC  2.  -1  KL 

Q  is  a  continuous  function  then  Q  (e  ,  C  ,  0)  -  0  .  A  similar 

2  3 

argument  being  valid  for  Q  and  Q  we  have  the  proof  cf  (5*25). 

The  constitutive  functions  may  also  be  expressed  in  terms  of 

-Ikl 

C  .  .  For  this  we  note  that  C  is  invertible.  Hence  we  may  write 
KL 


(3.26)  --  *o(e  ,  Ckl) 


(3.27'  «K  *  ¥e  '  °KL  '  V 


(5-28)  x“K  x“L 


Equation  (3*21)  for  the  stress  tensor  becomes  in  this  case: 
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'3.29) 


,  a  .fi 

iaQ  x  K  >  j 


The  argument  for  material  symmetry  restrictions  follows  the 
classical  lines  for  various  crystal  classes.  In  particular,  for  an 
isotropic  material  with  a  center  of  symmetry  is  a  scalar 
invariant  under  the  full  orthogonal  group  of  transformations,  {£}  , 


qkl 


The  free  energy  is  thus  a  function  cf  the  three  invariants 

-1 

I  =  tr  £  =  tr  £ 


(3.30)  I, 


"■*-  p  O 

tr(C)  =  tr  £ 


-1 


lx  s  tr(£)5  =  tr  £ 


For  an  isotropic  material  ifQ  can  hr  c^usidered  a  function  of 
8  and  Cqp. 


♦o  =  +o<6  '  V 


One  can  also  show  that  (3.21)  is  equivalent  to 


(5'51)  =  -2TJ0  ^  °7P 


From  the  reduced  form  of  the  Clayley -Hamilton  theorem  (A.4l) 


we  have 


.  .2  .  .  (I15+2IJ-3I2  Il)  „ 

zi  £  — 5 —  £ + - k - a 


The  stress  tensor  is  finally 

op 


(3.52)  t. 


o  2 


-“o(Ix  +  2I3  '  3I2  h)  5T  sce  -  no[25i7  -Xh  -htt 
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■no[4  *r2  *  6  xi  3^  cor  crP 


For  isotropic  materials,  following  a  similar  argument  to  that  of 
nonrelativistic  continuum  mechanics  [33],  it  can  he  shown  that 


(3.33)  qp  «  (Kx  Bap  ^c3s  +  K3  c“r  e^)  ( +  •  Dua) 


where  and  are  scalar  functions  of  the  invariants 


Z1  *  X2  '  I3 


(3.34) 


*  *6 


*08*  „  *  a  Qy  * 

ea  0  V  e5  ‘  ea  c  p  c  er 


The  inequality  (3*24)  becomes  in  the  isotropic  cac.e 


(3.35)  ei  +  *2  03  +  K3  03  1  0 


The  only  other  theory  of  this  order  of  generality  is  that  of 
Bresaan  [24].  The  stress -strain  lav  which  was  deduced  in  this  section 
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froa  the  second  law  of  thermodynamics  was  arrived  at  by  Bressan  by 
assuming  that  a  stress  function  exista.  The  heat  conduction  law  (5.33) 
is  a  generalization  of  the  one  proposed  by  Eckart  (17).  Bressan, 
as  does  Pham  Man  Quan  [25i,  attempts  to  generalize  Fourier's  law  of 
heat  conduction  without  any  regard  for  the  second  law.  It  seems  that 
it  is  impossible  to  satisfy  the  second  law  (as  stated  by  (2.52)) 
without  eliminating  the  heat  flux  if  one  assumes  their  foimj  of  the 
heat  conduction  law. 


Si 


framework  the  dependence  of  the  constitutive  equation*  on  i» 

a  constitutive  assumption  defining  a  class  of  materials. 

The  procedure  followed  here  is  identical  to  that  used  for 
thermoelastic  solids.  The  second  law  of  thermodynamics  for  non¬ 
polar  materials  reduces  to: 


n  cfy 
o  /  o 

?■ 


q  ^  n  blf  u  ndijr  * 


(5.37) 


P  6 


6 


2  c®  v  * 

+  no  ^  8  ]  % 


6 


>  o 


dd 


Q0 


where  we  used  (2.l4)  and 


At  any  event  £  the  following  quantities 


(3.38) 


(1,J  =  1,2,3) 


can  be  varied  independently  in  the  sense  described  by  (3-10).  By 
repeated  use  of  the  argument  employed  to  deduce  (3.14)  from  (3.13) 
one  finds  that  a  necessary  and  sufficient  condition  that  the  second 
law  (3*37)  is  satisfied  is  that1 

(5-59)  n00  =  5^ 


1Since  all  the  components  of  the  tensors  in  the  set  (3.36)  are  not  in¬ 
dependent.  in  deducing  (3.40)  to  (3*42)  one  should  consider  V  as  a 

*  *  u 
function  of  B  ,  nQ  ,  d^  ,  0^  ,  . 


0 


(3.40) 

S 

(3.4l) 

“T" 

(3  42) 

(3.43) 

6 

where 

(3  M) 

=  0 


0 


s 


0 


is  the  dissipative  stress  tensor. 

By  using  the  results  given  in  the  appendix  on  matrix  functions 
of  Lorentz  invariant  functions  and  the  conclusions  (3-39  to  (3.^) 
for  the  satisfaction  of  the  second  law,  one  deduces  that: 

The  constitutive  equations  of  a  thermo viscous  fluid  satisfy 
the  second  law  or  z he rmodynami c s  and  the  requirement  of  Loire ntz 
invariance  if  and  only  .If 

i.  The  free  energy  depends  only  on  the  temperature  and 
particle  number. 


(5.45)  =  *0(6  ,  nQ) 
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(3.46) 


11.  The  «ntrppy  and  stress  tensor  a r*j  63  tc  rained  from; 


oo 


0 

3<r 


2  C^' 

t  s  -  n  S  +  t 

*  o  on  *  D* 

o 


ill.  frhe  energy  flux  q  and  the  dissipative  part  of  the  stress 

1 


tensor  pi  have  the  form : ' 


*2 


#  * 


.*  *  *  * 


1)1  :  +  +  ^  +  \  k  Q  £  +  ^t£  ®4£  +  4£®£] 


.###*. 


,*  #p  *  #2*  * 

+  X6(£®<T£  +  4  £®£]  +  yjfcxi -4*£] 


(3.47) 


.*  -»o  #o  * . 


+  ^8tjg  x  ic  -  s?  x  £]  +  yid  x  l2)  -  ihl  x  I)] 


and 


+  \0d®dx4£>  ►  (|x4|)8il 
+  *ud  ®(£  x  42  £)  *  <£x!2£)®£] 


(3.M) 


*  *  *  *0  *  *  *  *  #*o# 

^  fl  +  n2  4  £  +  ^  4  £+y^4£+yx4  £ 

.*  *  .*2  *v 
+  ^(4£)x(4  £) 


where  and  are  functions  of  the  Invariants 


Id  .  tr  4  ,  IId  >  tr  42  ,  IIId  .  tr  45  ,  ^  .  £  •  £  , 
d  *  *  *  d  *  *2  *  d  *  ,*  *v  ,*2  *v 

ec  .  £  4  £  ,  e5  »  £  4  £  ,  \  -  £  1  (ft  £)x(4  £) 


where  pi  and  q  must  satisfy  the  Inequality  (3.43 ). 


have  not  assumed  that  the  fluid  possesses  a  center  of  symmetry.  In 
the  latter  case  further  reductions  are  possible.  *‘or  this  use  conclusions 
made  in  the  appendix. 


CHAPTER  IV 


EI£CTROMAGHETIC  THEORY 

!Hie  tias  1  c  laws  of  electromagnetic  theory  are  the  conservation 
of  charge,  the  conservation  of  magnetic  flux,  and  Ampere 'b  and 
Guass'  Lave.  The  relativistic  formulation  of  these  lavs  is  well 
known  (cf.  Post  [54],  Truesdell  and  Toupin  [8]  and  Miller  [l4J). 


Conservation  of  Charge 

To  formulate  the  lav  of  conservation  of  charge  assign  to  every 
three  dimensional  suhspace  in  space -time  a  scalar  function 
of  the  form; 

(4.1)  Q[s5J  «  J  oa  ds^ 

B5 

vhe_*  o  ia  called  charge  -current  vector.  The  lav  of  conservation 
of  charge  states  that  Q[s^]  vanishes  for  every  three  dimensional 
circuit. 

(4.2)  j  o“  =  0 


By  familiar  arguments  of  the  Green-Gauss  theorem  this  leads  to  the 
following  differential  equation  and  Jump  conditions. 
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(4.5) 


a 

0  .a  B  0 


,  (<ja]  Z  =  0  across  Z(x^)  •  0 


a 


It  is  useful  to  resolve  the  vector  a  into  ft  sp&ce-like  vector 
and  a  tin* -like  vector,  i.e., 


a 

a  £  -  o  u 


a 


(*»  /  8^oa  ,  -  o 


o  z  n  a 
o  o 


so  that 


(4,5)  a*  «  n  a  ua  +  ja 

n  O 


Equation  (4.5)^  with  the  use  of  (2.14- )  becomes 


(4.6)  nQ  Dao  +  Q  «  0 

The  physical  interpretations  of  a  and  are: 

cr  is  the  charge  density  in  the  instantaneous  local  rest  frame, 
a  r  i  i 

J  reduces  to  IJ  ,  OJ  in  the  instantaneous  local  rest  frame  where 
J1  is  the  conduction  current. 


Conservation  of  Magnetic  Flux 


The  conservation  of  magnetic  flux  is  obtained  by  assigning  to 
every  two  dimensional  subspace  in  space -time  a  scalar  quantity 
9  [SgJ  ,  called  the  magnetic  flux.1 


(4.7) 


$ [s_ ]  5  i  dsJ*5  :  e  _  =  -4' 

1  2J  2  J  op  2  *  op  ( 


The  lav  of  conservation  of  magnetic  flux  is  the  statement  that  the 
scalar  #[s gi  vanishes  for  every  two  dimensional  circuit. 


(4.8)  j 


•ce4^0*  =  0 


Ey  use  of  Stokes1  theorem,  this  leads  to  the  well  known  equations 


(4.9) 


*  «  0 

7&,P 


^Qp7b  ^  j  j  *  o  across  ^(x11)  »  0 
—  70-  ,p 


The  tensor  4>  is  physically 
op 


(4  10}  ^  =  (dual  £  ,  JS] 


where  JJ  i*  the  density  of  magnetic  flux  and  g  is  the  electric  field. 


1Ths  surface  element  of  the  surface  with  parameterization  u^,u^ 

.......  op  „  ! 


is  defined  by  d*2  “  2  au  dul  du2  * 
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That  la: 


(4.11) 


Maxwell* ■  equations  (4.9)  are  often  encountered  in  the  literature  in 
the  form 


(4.12) 


♦  _  +  ♦  +  ♦ 

op, 7  $7, cc  7a,  e 


o 


This  io  obtained  by  multiplying  (4.9)1  by  and  using  equation 

(A. 55).  Sometime b  it  it  convenient  to  use  the  dual  of  ♦  ,  ♦ 


(4.15) 


=  -  i  c**  ♦  ^ 

2  76 


In  term*  of  £  and  £  , 


A 

♦ 


is : 


c  [dual  £  ,  £3 


The  con*ervation  of  magnetic  flux  (4.9)  in  term*  of 


A 

♦ 


reduce*  to: 


(4.14)  S'*  «  0  ,  [♦*]  z=o 

>P  *■"  "  }P 

It  i*  poaeible  to  d^compoae  ♦  (cf.  Jailer  [l4 ] )  into 


6B 


(^15)  f<*  *  £  0  ua  -  2  a  ue  +  \* 


where 


£a  *  >  ^a"“  =  3  » 


(4.16)  *  ,5 

E  Sa  Sp  %&  ; 


*  8 

;  ^  up  »  0 


The  spatial  tensor  is  conveniently  written  as 


5  <W*  ®  7  u&  »  Pa  u„  :  o 


a£  “  Q876 


(4.17) 


©a  1  Q07&  A  1 

"  2€  9  2 


1  ***  * 

2  V  ^ 


The  magnetic  flux  tensor  $  _  becomes: 

ap 


(U.18)  v  .  ep«a-  $aup  +  caprB  e7  u 


'  '  uu 

ar'd  the  dual  of  £  ,  ♦  ,  can  he  decomposed  into: 


(4.19)  V*  =  ©auP  -  © p  ua  ♦  «<**£  r  u, 


Ihe  two  (^composition  (4.l8)  and  (4.19)  are  useful  for  the  formulation 
of  constitutive  equations  for  the  electromagnetic  quantities  and  for 
expressing  the  interaction  of  electromagnetic  fields  with  matter.  In 
terms  of  E  and  £  the  four  vector  £  a  and  ©a 


are 
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(4.20) 

fl  +  zx£  X  •  £  ] 

!&a  .  r«  - 1*  £ 

,  £JJ] 

*Jl  -  V2  Vl-v2 

Vl-v2 

Upon  substituting  (4.19)  into  (4.14^  and  multiplying  the  results  by 
uQ  and  S7a  and  using  (A. 26)  we  find  the  following  equations  which 
are  equivalent  to  (4.14 )^. 

(».a)  s\  8’,,  -  •*"  $  „  ■  0 

<»•■)  !*'  S,J#  •  ?*’  &  »,  *  »SB  *  8“  -  s'  -  » 

It  should  be  observed  that  only  three  of  the  fair  equations  of  (4,22) 
are  independent. 

Aripere*s  and  Gauss*  Lavs 

Andre's  and  Gauss*  lavs  are  combined  into  one  invariant  lav  by 
assigning  to  every  tvo  dimensional  subspace  s^  a  scalar  invariant 
r(s2] 


(4.25) 


r(s2] 


•/ 


<3$  A 
G  ds, 


& 


200 


^  ;  **2c$  S  2  *0^76  **2 
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Ampere *s  and  Gauss*  laws  state  that  for  every  closed  two  dimensional 
circuit  enclosing  a  three  dimensional  subspace  s^ 


■jiMnwwar 


TO 

(4.24)  /  0°^  d.^  --  f  oa  ^ 

The  differential  form  of  (4.24)  and  the  Jump  condition#  are: 

(4.25)  0°^  a  s  °a  t  £g°^]  Z  ■  0  across  X(x^)  s  0 

,P  ~  -  »p 

os 

The  tensor  G  is  the  electric  displacement -magnetic  field  intensity 
tensor : 

(4.26)  G°^  =  (dual  g  ,  -g] 

where  £  is  the  electric  displacement  and  g  is  the  magnetic  field 
intensity.  Explicitly  (4.26)  gives 

(4.27)  01J  =  e1*  H,  .  G4*  -  -Oi4  =  D1  ,  G44  »  0 

Similar  to  (4.15)  ve  decompose  G^  into 

(4.28)  G0^  .  t)  P  ua  -  Va  up  +  s0^76#  7  ^ 
where 

£>p  =  0^^  ,  £>Pup=0  ,  ^pUP  >  0 

(4.29)  G06  =  Sar8P8G76=  c0676  H  7  u6 

'  |  <*»  f  *  ■  i  97  »' 
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In  tern*  of  the  electromagnetic  field*  £  and  £  ,  and  are; 


(*.30) 


.  r£  +  X  X  I 

/l  -  V2 

x  •  £ 

'/l-V2 

X  13 

„  rfi-zxfi 
Vi  -  y2 

x  •  S 

’  3Z* 

By  substituting  (**.28)  into  (**.25 )1  in  the  same  vay  as  done  in  obtaining 
(**.2l)  and  (4.22)  ve  get 


(4.31) 

o’,,  * 

(4.32) 

Only  three  of  the  four  equations  (4.52)  are  independent. 


CHAPTER  V 


i  «■  minm— ilifWf' 
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ELECTROMAGNETIC  INTERACTIONS  WITH  PONDERABLE  MATTER 

The  presence  of  matter  in  an  electromagnetic  field  has  been  the 
object  of  researches  since  the  beginning  of  electromagnetic  theory. 

There  are  several  approaches  to  this  problem.  The  one  that  is  adopted 
in  this  chapter  was  originated  by  Lorentz  [553  in  order  to  derive  the 
electromagnetic  field  equations  for  ponderable  matter.  According  to 
this  approach,  the  interaction  between  matter  and  the  electric  fields 
can  be  deduced  from  a  microscopic  model  to  within  the  order  of  approxi¬ 
mation  desired,  The  form  of  the  interactions  so  deduced  can  then  be 
taken  as  the  starting  point  for  the  development  of  a  continuum  theory. 

The  force  on  a  material  body  can  also  be  determined  by  assuming  an 
effective  current  distribution  and  postulating  that  the  body  force  is 
the  Lorentz  force  on  this  distribution  of  charge -current .  By  redefining 
the  energy-momentum  tensor  one  can  show  that  these  two  approaches  are 
equivalent.  A  third  approach  to  this  problem  is  to  attempt  directly 
to  write  down  an  energy -momentum  for  the  material  body  and  the  electric 
field.  Such  attempts  are  usually  guided  by  one  of  the  above  approaches. 
There  is  still  no  widely  accepted  form  of  the  energy -momentum  tensor. 

The  point  of  view  adopted  in  this  chapter  is  that  matter  is  acted 
upon  by  forces  of  various  types,  one  of  which  is  due  to  the  electromagnetic 
Interaction  with  the  molecular  and  atomic  structure  of  matter.  By 
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applying  the  arguments  of  Dixon  and  Eringen  [6],  it  is  physically 
reasonable  to  assume  that  electromagnetic  interaction  with  matter 

(neglecting  electric  quadrupolea  and  higher  terms)  is  due  to  a  body 

^  1 
force 


(5.1)  f/  =  ^ 

and  a  body  couple 

(5.2)  L^V  »  -  irj*1 

06 

where  T r  is  defined  by 

(5.5)  «  ♦*  -  ^ 

06 

The  polarization  tensor  tt  is  a  skew-symmetric  tensor  of  the  form 
(5.4)  irQp  =  [dual  g  +  dual  (X  x  £)  ,  £] 

where  £  is  the  polarization  vector  and  JJ  is  the  magnetization  vector. 

This  definition  of  tt0^  is  that  of  Lorentz.  It  is  well  known  that  it 
lacks  a  symmetry  in  the  transformation  of  the  magnetization  vector  and 

1In  body  force  T±  of  [6,  eq.  5.21]  we  drop  the  quadrupole  moment  tensor 
q1^  and  take  £  *  £  -  d(p  X  k)/dt  where  £  is  that  used  by  Dixon  and 
Eringen.  This  rate  term  will  appear  later  in  the  momentum,  cf.  eq,  (5.18) 


below. 
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The  magnetic  term  in  the  electromagnetic  body  force  (5*1)  if 
based  on  the  Amperlan  current  model.  Some  researchers  prefer  a 
magnetic  dipole  model  (cf.  Fano,  Chu  and  Adler  [26]  and  Fenfield 
and  Haus  [36]).  If  one  usee  this  model  then  instead  of  (5.1)  one  vculd 


assums 

f/  *  U°  Up)  ^  ua)i^p  +  a“  ^ 

where  and  <i6a  are  defined  by  (5*8)  and  (5 .9)  respectively. 

Since  the  existence  of  magnetic  poles  is  doubted  we  prefer  the  Amperlan 
model.  It  should  be  noted  that  both  (5*1)  and  the  form  listed  in  this 
footucte  for  the  magnetic  dipole  model  lead  to  the  same  form  of  the 
first  lav  of  thermodynamics  (5.22).  Thus  the  conclusions  of  Chapter  VI 
remain  valid  for  the  magnetic  dipole  model. 


the  polarization  vector.  This  has  no  effect  on  electromagnetic 
phenomena. 

The  balance  equations  (2.25)  and  (2.^0)  are  now 


(5.5) 

** 

JP 

(5.6) 

TUv] 

>> 

♦ 

:± 

“m® 

II 

+  0°  *1*  +  f“ 


It  is  useful  to  decompose  t3®  into 


(5.7)  tP  *  f  f  6  ua  +  5* 


where 


*-  TT^U 


a 


?*  =  sa  s\  /6 

7  & 


(5.8) 


f 


75 


One  can  also  define 

(3.9)  Ma  --  |  ‘op^  f1  u8  ■ 

Brno 

(5.10)  £**  =  s^75  ^  r  Ug 

The  polarization  tensor  tP  reduce#  to 

(5.11)  T«  = 

It  should  be  noted  that 

(5.12)  9#  Up  =  0  ,  ,4(,P  Up  =  0  ,  Up  »  0 

Using  (4.l8)  and  (4.28)  ve  see  that  (5*5)  is  equivalent  to 

(5.15)  Va  -  fa  +  Pa  ,  *a  =  6a- 

The  spatial  part  of  the  four  vector  “P  ^  reduces  to  the  polari¬ 
zation  vector  and  the  spatial  part  of  tM-  ^  to  the  magnetization 
vector  in  the  local  instantaneous  rest  frame. 


7 6 


nomww 


This  interaction  model  la  a  physically  reasonable  description 
of  a  charged,  conducting,  polarizable,  magnetizable  material.  Hie 
following  observations  allow  the  construction  of  various  subclasses 
of  materials. 

a)  If  aa  ufl*o  ,  the  material  is  charge -free  in  the 
local  instantaneous  rest  frame. 

b)  If  S°.  a*3  *  ja  »  0  ,  the  material  is  a  non-conductor 

P 

in  the  local  Instantaneous  rest  frame. 

c)  If  7?^  s  0  (*Xa  *  0)  ,  the  material  is  unmagnetized 
in  the  local  instantaneous  rest  frame. 

a 

d)  If  a  o  t  the  material  is  unpolarized  in  tho 

local  Instantaneous  rest  frame. 

From  (5.7)  and  (4.18)  one  can  show  that 

(5.14) 

-•  Atf^ua  +  sa  .  &*oi 

where  (4.l6),  (5-12)  and  (A. 56)^  of  the  appendix  were  used.  Therefore 
the  right  hand  side  of  (5*2)  becomes: 

(5.15)  -  Tft[a  .  rlaS  ^  u3'1  +  *  uo]  +  ^ 

p  p  p 

From  (2.25)  the  left  hand  side  of  (5*2)  is 
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(5.16)  T[<*]  -  u[V3  +  pCau0]  -t1*1 


Equation  (5*6),  using  (5*15)  and  (5»l6)>  lead*  to 


(5.17)  P[a  Ue]  -  qta  -  t[ae3  ~  *>  up3  +  T  ^  u?3 


-  ^a<5&1  -  ^[aep] 


By  taking  the  projection  of  (5.17)  along  uQ  ,  one  obtain* 


(5.18)  pa  .  q“  +  r7  J7“ *  S7Pa 


Applying  the  projector  8  to  (5. 17),  one  deduces  that 


(5.19)  t[apl  -  f[a  +  Jl}a  ©Pl 


Therefore  a  necessary  condition  that  the  balance  of  angular  momentum 
(5.2)  is  satisfied  is  that  (5.18)  and  (5*19)  are  valid.  It  is  easily 
shown  that  (5«l8)  and  (5*19)  are  also  sufficient  for  (5*2). 

To  obtain  the  first  law  of  thermodynamics,  we  substitute  tj1  +  f^ 

(where  is  given  by  (5*l) )  into  (2.51).  Hence 

n  Dc  +  +  pa  Du  -  t°^  u  -  tt®  R  u  -  c*  ♦  u  -  f^  u 

o  ,P  a  op  a,p  p.  a  »  ^ 


(5.20) 


-  -  ??»&a  -  s>a»pan»p  - 
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Where  (4.4),  (4,l6),  end  the  identity 


(5.21)  f*  ♦^p  -  vea  *  $a  ¥*  DUp 

have  been  engjloyed.  To  prove  (5.21),  one  uses  (4.15)  and  (5.7)*  With 
the  help  of  (4.l6)  and  (5- 12), the  first  tvo  terms  on  the  right  hand 
side  fall  out  almost  Immediately;  the  remaining  tvo  terms  are  arrived 
at  by  using  (A. 56)  of  the  appendix  and  (4.22). 

Substituting  (5*18)  and  (5.19)  into  (5.20)  ve  obtain 


(5*22) 


■v-  *  «s,s  * “j  - «<®1  **  -  iP,I<V1  »«,„ 
•  a‘“  e el  vs  „ »  »e0 .  fs , 


Equation  (5.22)  is  the  first  lav  of  thermodynamics  for  an  electro¬ 
mechanical  material.  It  is  often  convenient  to  revrite  (5.22)  as: 


no  ft  ♦  /p  ♦  DUp  -  [*<*>  -  $  P)  -  UL{a  8  P)] 

(5.25)  +  <pa  ID  S  a  -  5  P  «a,p)  +  Jf  l»S  a  -  *6°  ua>p] 

-J 


An  equivalent  third  form  of  (5-22)  is 
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n  D[c  +  fp.  ♦  -  ♦  i  Du. 

0  nQ  nQ  >P  P 

(5.24)  -  it(ae)  +  r(a  8  e)  ♦  M.(ae> 6i  +  s°p  +  r7<5  s01^]! 

7  7  ~ 

-  8a  [DSPa  -  F*  %g]  -  »“  &•*«  -  ua>jJ]  -J "g  a  -  -f“u 


Each  of  these  form*  may  he  found  useful  in  various  situations.  For 
example  (5.22)  is  convenient  if  one  wishes  to  describe  a  fluid.  The 
form  (5.23)  will  he  used  for  a  solid  when  the  independent  variables 
are  &  a  and  #3  a  .  If  the  independent  variables  are  9°  ^  and 
ciC*  >  then  (5.24)  is  the  appropriate  fora  of  the  first  lav. 
According  to  the  entropy  inequality  (2.51) 

P  f *  •  \i 

(5-25)  nQ  D^o  +  (a_)j0  +  __Ji  >  0 


Employing  (5.22),  (5*25)  becomes 


(5.26) 


“o<D,1oo  -  ?  E‘  >  -  £  V 


t(°e)  .  r[a<g pl  . 

■  d  +  mi.  i.  1  1  ■  fa) 

•3  ap  e  ap 


M[a&*]  *  ^<?a  jfaea  f&a  ^ 

9  wap  "  9  ~  9  "  +  6  - 


where  is  defined  by  (5.2). 

Using  (5.23)  in  (5.25),  an  equivalent  expression  to  { 5 « 26 )  useful 


fo.  solids  is 


8o 


(5.27) 


Another  form  of  the  second  law  is  possible  if  one  uses  (5*24). 

Using  (4. 25)  and  (5.5)  to  eliminate  aa  from  (5.1)  *nd  employing 
(4.12),  it  is  possible  to  express  the  body  force  as  the  divergence 
of  a  tensor 


(5.28) 


V6  = 


n  J  1  .OT  t  jiP 
7°  ‘IT*  *7Vr 


In  the  derivation  of  the  Jump  conditions  (2.25)  it  was  implicitly 
assumed  that  body  force  is  continuous  across  the  singular  surface  2  . 
This  is  not  necessarily  the  case  for  (5-1)  since  the  electromagnetic 
fields  may  suffer  a  discontinuity.  The  forms  of  the  balance  of  energy 
momentum  and  its  Jump  conditions,  taking  into  account  this  possibility, 
a  ie 


(5.29)  (T^  +  T  ^l3)  .  -  f^  ,  [T^  +  T  1 2—0  across  2(xM)  =  0 

e  ,p  ~  e  —  ,p 

The  balance  of  i.  ment  of  energy -momentum  (5*6)  becomes: 

(5.50)  +  T  ^  = 


0 
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Thus  it  is  possible  to  /rite  a  "total"  energy-momentum  tensor  for  the 
interaction  model  presented  in  this  chapter.  This  tensor  is  symmetric 
as  a  consequence  of  the  principle  of  balance  of  moment  of  energy- 
momentum.  Whether  it  is  possible  to  obtain  such  a  tensor  for  any 
interaction  model  is  an  open  question.  The  important  point  to  remember 
is  the  Jump  condition  (5*29)2*  Whether  the  Maxwell  stresses  are  "real" 
stresses  is  a  much  discussed  point.  It  seems  unlikely  that  they  are 
(see  Dixon  and  Eringen  [6]  for  a  discussion  of  this  point).  The 
introduction  of  the  interaction  energy-momentum  tensor  is  a  mathematical 
convenience.  The  system  consisting  of  the  body  forces  (5*l)  and  the 
body  couple  (5*2)  is  equipollent  to  a  system  of  "surface  tractions" 
given  by  (5.28). 
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CHAPTER  VI 

CONSTITUTIVE  THEORY  OF  DEFORMABLE  ELECTROMAGNETIC  MATERIALS 

A  constitutive  theory  of  deformable  electromagnetic  materials 
can  be  formulated  following  the  principles  enunciated  in  Chapter  III 
In  this  chapter  relativistic  theories  of  polarizable,  magnetizable, 
conducting  solids  and  fluids  are  presented.  The  requirement  that 
the  constitutive  theory  be  thermodynamically  admissible  leads  to  a 
considerable  reduction  in  the  form  of  the  constitutive  equations. 

The  following  theories  are  sufficiently  general  to  include  thermal, 
electrical,  and  mechanical  effects  with  restricted  spatial  and 
temporal  variations.  Thus  such  effects  as  gyrotroplc  phenomena, 
optical  activity  and  heredity  are  excluded. 

Electromagnetic  Solid 

For  an  elastic,  magnetized  dielectric  which  is  also  a  conductor 
an  appropriable  set  of  independent  variables  is : 

(6.1)  e  ,  xpK  ,  /  ,  <fK  ,  SK 

JC  K  v 

where  v  ,  S  and  o  arts  defined  as: 
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(6.2) 


<gK  -  *gn  (xi/xK) 

>P 


The  sgn  (X  /^K)  ,  which  signifies  the  sign  of  the  Jacob ian,  is  introduced 
so  that  ©  ^  remains  as  an  axial  vector.  The  dependent  variables  are 


(6.3) 


t(oe) 


.a 


with  p°  and  t^°^  determined  from  (5*18)  and  (5.19)  • 

One  can  now  proceed  in  a  manner  identical  to  that  used  in  Chapter  III 
for  mechanical  materials  s  write  the  entropy  production  for  this  set 
of  constitutive  equations  and  find  the  necessary  and  sufficient  condi¬ 
tions  for  it  to  be  non -negative.  To  this  end  we  use  the  chain  rule 
of  differentiation  in  calculating  D«  and  Dnoo  and 

D<gK  -  X^p(Dg  P  -  Qa  uP>a) 

D«K  .  /.(DC*3  -  ©auP  J 

tP 

Recalling  the  identity  (1.55),  the  second  law  ( 5»27)  becomes 


(6.4) 


n  dt  d* 

— -  (^-2.  +  tj  )  De  -  i-  (a  — 2- 
S  oe  oo  8  o 

.[t(o3)  _*>(«^  P)  -  Ji(a  8^1  7^  X*  p)  DX* 

n°  i/  ^  4-  \  r,/PK 

--  JJtM  -  5<  PK  ‘  %  5T>  '« 


Sc 


«F 


>  0 


where  K  *  JK  and  ^ 


are  defined  as 


r. 


ji 


=  ?°  a 


a  x  K 


=  Jl  xa„  e  'x 


(6.5) 


'K 


a  K 


Jax0K 


gnr/jK) 


Sc  '  ^ 


a 


K 


and  +q  is  the  free  energy  tQ  *  «  -  8  *1Q0  .  At  any  point,  £ 

space -time  the  following  quantities  can  be  assigned  arbitrarily: 


(6.6) 


e  ,  ,  3K  ,  8K  ,  De  ,  dxpk  , 


»eK  ,  d SK  ,  d8 
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,  in 


Thus  by  the  argument  leading  from  (5*15)  to  (3.14),  a  necessary  and 
sufficient  condition  that  the  entropy  production  (6.4)  is  non -negative 
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la  that  the  following  idantitiea  hold: 


(6.7)  n 


00 


*o 

3iT 


(6.8) 

(6.9) 

(6.10) 

(6.11) 

(6.12) 

(6.15) 

(6.14) 


t(oe)  .  8r(a  P  )  +  0)  +  ^(a  3  P) 

<*  K 


K  0  d£K 

^  ='n°SK 


■  0 


*o 

a7 

V  -  o 

s7(axP)  fjo  .  0 

^  Av7 

K 

0  flK 

V  *  0 


Aa  in  Chapter  III,  In  order  that  t>*  invariant  under  the  proper 
Lorentz  group,  fQ  must  he  a  function  of  the  form: 


(6.15)  ♦„  ■  t0(»  ,  Cffl  ,  £K,  ©K) 


"VW*1 
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8ince 


*o  .  *0  *0  0  *o  0 

^  >5  ^  7*^» 

equation  (6.12)  and  (6.13)  are  satisfied  Identically  and  (6.8)  becomes 


(6.16)  t(ap)  -  2n0^K^L  J2-  +  P(a^P)+  ^(“»P) 

KL 


In  a  similar  way  and  must  have  the  forms 


J.  *  Me  ,  P  ,  &K  ,  &K  ,  cn) 


(6.17) 


Sc  =  V®  >  ^  ’  8  K  »  ®K  ’  Ckl) 


We  have  thus  arrived  at  the  following  important  conclusion: 

The  most  general  form  of  the  constitutive  equations  for  an  aniso¬ 
tropic  elastic  dielectric  with  heat  conduction  satisfying  the  require¬ 
ment  of  Lorentz  invariance,  non -negative  entropy  production,  and  the 
balance  of  moment  of  e nergy-mome nturn  is  the  following: 
i .  The  free  energy  assumes  the  form 

(6.10)  t0  ■  t0(* ,  C  ,  8  K  ,  «K) 


06 

ii .  The_  entropy  H00  ,  the  stress  tensor  t  ,  the 
polarization  vector  a  j  and  the  magnetization 
vector  M  a  are  determined  from  the  free  energy  by 
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(6-l9)  \o  =  -  ST 

(6.20)  t<*  *  2nfl  *V@L  9*  e. 15  ♦  ^°sp 

KX 

(6.21)  *  -nQ  -p, 

(6.22)  Ma  -  -n0  f  a  *80  (VyK) 

vhere  (5*19)  employed, 

ill.  The  nonmechanical  momentum  £  1b  determined  through  (5.18) 

(6.23)  p  5  1  ?  x  ®  -  §  x 


iv.  The  conduction  current  Ja  and  the  heat  flov  vector 


are  determined  from  the  equations 


(6.24)  Ja  -  x  a  JK 

(6.25)  ^  -  X1^ 


vhere 


(6.26)  jk  -  JK(e  ,eK,gK,  gs  ,  ctt) 

(6.27)  Qjj-  -  ^(e  ,  ,  S  K  ,  6  K  ,  c^) 
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v.  The  conduction  current  and  the  heat  flow 

. . ■"  ■!  II  111  ■■■■!  «»  '  ft, 

vector  must  satisfy  the  inequality 

_  K  VK 

(6.28)  JK<?  -  -=5-  >  0 


which  in  view  of  the  continuity  of  JK  and 
implies 


(6.29)  JK(e  ,  o  ,  o  ,  gK  ,  cKL)  .  o 
(6.%)  ^(e  ,o,o,  gJK  ,  c1®1)  -  o 


Further  reduction  in  the  form  of  the  free  energy  >  the  heat  flow 
vector  ,  «.ad  the  conduction  current  can  be  obtained  if  one 
knows  the  crystal  class  of  the  material.  In  particular,  for  an 
isotropic  material  with  a  center  of  symmetry,  a  minimal  integrity 
basis  for  the  symmetric  tensor  <;  ,  the  vector  Q  ,  and  the 
axial  vector  ^3  K  (see  Smith  [33))  is: 


C*  I 
C  K  '  h 


c™'L  ,  I,  ->£_> 


MK 


MN  K 


(6.31) 


0  *  >  E! 


0>K~1  ff  L  *K;1  -in  a  L 

CKL  ®  ’  E2  &  CKM  l'  L  & 


(Eq.  continued  next  page) 

It  is  convenient  to  consider  yQ  as  a  function  of  instead  of 

C-_  .  The  form  of  (6.20)  corresponding  to  (3.21)  is  easy  to  deduce. 

KL 
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(Eq.  6.31)  continued) 


B0  -  #k6k  ,  Bx .  eK  c^1 ,  b2 . 6*  -4  i\eL 

bj  ■  c1^ 

\  -  (8KeV ,  b5  .  (gRe„)  U3K  4  gL) 

b6  -  (£H«R)  (©K  ckm  4s1) 

•  SKLh£K 
B9  ' 

B10  “  ^  €KIK  ^  °LI^  ® 

B11  -  <68©s)  C!CU|SK  ‘c1IR«R 


It  can  be  «ihovn  that  the  set  (6.31)  can  also  be  expressed  in  the 

fro*: 


tr  £ 


tr  £ 


tr  £** 


(6.32)  EQ 


£ 


2 


€  £  S  ,  \  -  $.£  I, 

o  (I^zyjv,)  (i12-x2) 

*  '  S - 5 - 2 - 


8  £  £  +  \  S<?  § 


(Equation  continued  next  page) 


(Equation  6.32  continued) 

B0  -§a«,  \  -ga’? 

b2  -«•  g^llY--1--  -  %*■  $+  \  ?*2® 

b5  -  |J|  »a(Sx  s2sb)  ,  ^  -  (ga©2  ,  Bj  -  (£a«)(2Jaa4) 

(1A21  .51 1  )  (1A1 ) 

b6-  (g  aS)  tS*«  — — V^~  -  ~V-g-g  &*♦ 
b7-  |j|  <a&)x«  -  g  ,  b8.  |j|  (£2&) 
b9  -  |J|  (a8)  x  a2  €  •  ?  ,  b10  -  |  jf  (g  c»)(a ©  x  5  • « 
Bu-(ga§)  (a2©  x?'S  |j| 

where  the  following  notation  has  been  used 

€aS  -  Sa^s\  (£<£)„«  =*«  p 

(6xg)a-  «aP7B®P<g7u6 


From  the  work  of  Vine  man  and  Pipkin  l3T)>  for  an  isotropic  material 
with  n  center  of  symmetry,  tQ  muat  be  a  function  of  the  invariante  (6.32) 


(6.33)  t0  «  t0(e 


/a -1,2,3  \ 

lb  -0,1,2  ) 

'C  —  0,1, a.., 11  * 
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From  (6.52)  after  some  algebraic  manipulations  ve  can  show  that  (6.20), 

(6.21),  (6.22)  assume  the  following  forms: 

ioT1£+T2£  +  T3£-2  +  TJ*£  £&  +  T5fe£*  £2  4  +  £2  S®  £4] 

+  g®  £  g  +  ^[£  §®  £2  ©+  SL  ©&  £  ©  ]+  *g[£  3?®(£x  £2  6) 

+  C?X£2g)®£«-  ©X£S)0£2S-£2e®  (®X£®' 

-  £g®  £(?X  £g)  -  £(§X  £g)  «  £  g]  +  ^[£80  £§  +  &§®£gJ 

(6.3*1-)  +  t1q[£2  £®£€  +  £22®£S+£§®£2S+£S«£26; 

-  ^[£g®  (5x§)  +  (gx?)®£g]  -  XgUgxS)®^ 

+  £4®  £(Sxg)+  £2<g®  (?x4)  +  £(«x£)®£<g]  -  ^[(fix  £<g>Sfe2£ 

+  £2  4  &  (?X  £  £)  +  £  <2  0  £(g  X  £  6  )  +  £(g  X  £  £)  ®  £  & 

-  £  £  ®(®  X  £2  e )  -  (fX£2g)®  £*]  -  X10[£'ga  (jgxfi) 

+  (£  X  g)  ®  £  g]  -  x11[£  S®  £U  xg)+  £(fi  x§)®£« 

+  (€x«)0£2§+C2S0(|x3)]  +  g 

ff—Xj^  -  X2££'X5£2^-X4S-X5£S-X6£2S  -  X^lfiGSxfc) 

(6.55)  -®X£§]  -  Xe[£2(gx<g)  -gX£2<gJ  -  ^[^(gx  £4) 

-  £(£x  £2  $)  ]  -  X10  S  X  £  g  -  XjjL  Sx  £2  © 

*1*!  5  -  £  S  -  £2  S  “  P^[£(g  X  £2  5)  +  (£29>  (£  2) 

-  £2(3x  £®]  -  ^ <z£  -  ^  £2 4  -  x  £^  -  xd  £x  £2  £ 

-  X9  (£©  X  (£24)  -  x10[£(£xfl)  -  £x  £fl] 

-  Xnf£2(^  X  §)  -  £x  c2  g] 


(6.56) 
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vher« 


(6.37) 


(L3+2I,-3T  I .)  3t  ftn  P  ftB 

"n0 - 1 -  37;  '  t2  -  -^o  31J  +  no{h  'V  37J 


3* 

-2n  ^-2n 


**-  **o 


-2n 


O  3^  -  “o  *1  3i^  ’  T4  ■  -^o  3I£  '  T3  “  '^o  3e^ 

at  at  at 

’  t7  *  -2no  sr  »  *8  “  ■  |j|no  3T 


at 


a* 


-“o«*»af  ,  T10.-no(€£S)^ 

(^21,-3121,)  *0  _  _  *„  .  „  2  T  ,  *o 

no - r - -  ai£  *  *2  ■  ar  ■  no(Ii  -I2)  ^ 

»>Wi> «, ,» . 


at. 


at 


35^  +  2noIl  31£  '  \  *  n0  6 


-  £  ~  3EJT 


at. 


at. 


b 


2x.0  «  £  +  n0(3£2«.)  3^  +  no  aj  U^tM^NMU 

(l13+2Ia-3I0I1 )  at 

+  ^— --  g  •  sj  +  n0  I  J|  «  .  U  §)  x  g  ^2_ 


+  nQ  !j|  (fi2#.)  x  gggjj 


at. 


at. 


at 


at. 


*6 

*8 


+  *7“a0  |j| 

no  |j|  3ET  ’  *9  ”  n0  |j|  ^,111D  =  -no  ,j|  3b^§£  6 


\L1 


atn  (i  3+2i,-}i„i. )  at 

-n0  |j|  §£§35-7  '  “l  “  no - r - 


o 


3+ 


O  3*- 

d  „  \  o 


^4 


3y  3+ 

3F  -  °o(Ii‘i'I2)  3^’^'2iio^+  “Vi  35“ 

3v_ 

no  35^  lJl 
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Also  using  the  results  of  Smith  [33],  we  can  show  that  the  conduction 
current  and  heat  flow  have  the  following  forms : 

J  -  “i  £  +  °2  £  s  +  °3  £2  S  +  %  &  +  °5  £  1  +  °6  £2  I +  aj  §  x  £  * 

®8  £<§*€)  +  ®9  £2(®x  €)  +  o10  gx  £  &+  ou  ?x  £2  ^  + 

S  x  £  +  o15  £(£(x  £)  +  ollt  £2(Sx  1)  +  Ojj  S  x  £  1  + 

(6.38)  „  *  D  9  3  » 

0^2  x  £  £  +  017[£(§X  £<£)  -  C2(gx  £%)]  +  Ol8[£(5x  £2  jg) 

-  £2(8  X  £  &)  ]  +  °19  2+  °20  £  ®  +  °2X  £2  °22  ®x  £  ® 

+  <*23  5x  £2  § 

q  -  Kj_  1  +  Kg  £  1  +  Xj  £2  1  +  ^  i  +  K5  £&  +  Kg  £2  £  +  ffx  1 

+  *8  £<ffx  1)  t  (C9  £2(2x  £)  +  *10  S?X  £  i  +  ltu  »  X  £2  l+KjjSxg 
39)  +  *14  X  <g)  +  *15  S  X  £  Sc  +  *#.*€♦  Kjjfoti&X  £2  Jg) 

-  £2(®X  £  £)]  +  Kjglzi&X  £2  £)  -  £2(gx  £  g)]  +  «19  3  ^*15£(®  x  £) 

+  *20  £  32  +  *2!  £2  6  +  <22  Sx  £  #+  *23^  *  £2  2 

where  a^g  and  k^,  k^q  aye  functions  of  the  invariants 

listed  in  (6.32)  and  the  following  invariants: 


* 


f 

I 

* 

\ 


£•£>£££>££  >  (i£©  GE?£  £)  > 

(1  £  S)  C&  •  &>>&  •  ®*  £  1  »  1  •  £2  1  ,  |  £(Cx  £2  £), 

(6.40)  &£©  (J5x  £*)  •  1  ,  (I  *  gx  £2  «)  C@£  1)  ,  S  •  1  , 

ffil,  ,  (#■€)  (i£®  ,  CP*  I)  (*£©  » 

Cgfig)  (§z&  ,  (la®  (®£2^)  ,  (i£®  («£2|)  . 

(I  •  ©  <»£2  I)  -  (I  •  ®  t?£2€)  ,  I  •  «X  &  ,  |  •  0X  £«  , 

(e  '  ?  x  £  I  »  I  '  ®X  £2  €  >  £  •  ®X  £2  |  ,  |  £(5 X  £2 «.) 

+  ie£(?x  £2  |)  ,  <i£®  @£(5x€)  +  (€£©S£(gx|)  , 

(I  •  *  x  £2  ©  («£€)  ♦  (§  ■  »X  £2«)  *?£  | 


and  ...,  have  the  form 

«19  “  (It®  “24  +  (§£®  “25  "  (®£2€)  “26  '  <§£2  |)®27 

“20  -  CS-  ®)  “28  +  (®'  I)  “29  +  (se£®  “50  +  (Ie®  “51 

(6,4!)  +  (S  £2<t)  “52  +  (P£2  I)  “35  +  ©£©x  |) 

+  ®  £  (0X«)  0J5  +  {<g  •  ?X  £2  §)  Ojg  +  (I  •  3x  £2S)  “57 

“21  *  Ufc©  “38  +  (-  £.©  “39  +  C?  •  ®)  “26  +  (£  •  S)  “27 

“go  “  -(I  £  5)  “34  -  (£  £  8)  '35 
“2J  -  (5  £  €)  “35  +  (*  £  I)  “57 


where  cr^  ,  ...  ,  are  functions  of  the  invariants  (6.5 2)  and 

(c.4o).  By  replacing  the  a  by  k  in  (6,4l)  one  obtains  the 

rs  r® 

functional  forms  of  k^,  ...,  . 


^•imw 


1 
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For  an  isotropic  material,  from  (6.32),  (6.35)  and  (6.34)  to 
(6.39)  the  set  of  dependent  quantities  <  ,  ij  >  %.  >  f?  >  *4? 

are  covariant,  tens  or  functions  of  £  ,  <2  ,  J0  and  j£  under  the 
full  Lorentz  group,  and  £  and  q  have  the  forms  (6.38)  and  (6.39) • 
(For  materials  with  no  center  of  symmetry,  they  are  corariant  tensor 
functions  under  the  proper  Lorentz  group.)  While  the  constitutive 
equations  for  isotropic  materials  are  deducible  from  the  general 
forms  (6.18)  to  (6.27),  in  many  cases  it  is  more  convenient  to  redsrlve 
these  equations  from  the  second  law  of  thermodynamics. 

The  second  law  (5*26)  for  an  isotropic  material  becomes  (using 
(1.55)  and  (1.65)) 1 


-  n  ♦  %  De  -  i  {t«*>  +  aptc  £  f>]  +  M  [a  e  Pi 


(6.kz) 


A  necessary  and  sufficient  condition  for  (6.42)  is 

"00  s  -3T 

(6M)  t<<*>  r  -  2n0  S°t  c’>  7 
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(6.45) 


(6.46) 


(6.47) 


(6.48) 


(6.49) 


s“t  [2n0  zrr  c°]7 '  ^[T  &  c]  -  **[T  &  a])  z 

7U  ' 


<pa  =  -t,  S“ 

V  0  7  dii^ 


a  ^0 


^  =  -V5  7  3F 


STT  V  UP  +  n0  JT  u6  £  a  +  n0  3T  UP  #<x  +  noTJ  s7a  z  0 


a  % 

i  fi  a  -  -r-  2  0 


The  requirement  that  ^  is  an  invariant  function  of  6  .  , 

o  Qp 

<?a  >  6a  and  uQ  under  Lorentz  transformations  gives : 

(6.50)  2^-ce]  y  +  t?-e Pl  +  +  =  0 


By  substituting  (6.46)  and  (6.47)  into  (6.50)  one  sees  that  (6.45) 
and  (6.48)  are  satisfied.  By  using  (5. 19)  in  conjunction  with  (6.50) 
it  can  be  established  that  a  necessary  and  sufficient  condition  that 
the  constitutive  equations  for  an  isotropic  material  are  Lorentz  in¬ 
variant  and  satisfy  the  second  lav  of  thermodynamics  is : 

i .  The  free  energy  t0  a  scalar  Invariant  function 
of  6  ,  £  ,  £  t  <8  i  U  ukfler  the  full  Lorentz  group, 
ii.  The  entropy  1qo  ,  the  stress  tensor  t  ,  the, 
polarization  vector  f  ,  and  the  magnetization 
vector  ^  are  determined  from  the  free  energy  by 


MwnMMn 
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St 


(6.31) 

*00 

'  mw 

(6.32) 

(6.33) 

T 

.  *0 

'  -no  *£ 

(6.5V) 

JC 

+*» 

.  *0 

“  "no  dg 

ill. 

The 

nonmechanical  p>  is  determined  hy  (6.23). 

iv. 

The 

conduction  current  J  and  heat  flow  vector 

<1  have  the  fora*  (6.38)  and  (6.39)  and.  satisfy 
the  inequality 


(6.55)  J  •  |  - 


a*  £ 

e 


>  c 


It  should  he  noted  that  (5.19)  i»  satisfied  identically  if  i.  and  ii, 
hold.  If  tQ  is  a  scalar  invariant  of  £  ,  <*  ,  i?>  J£  under  the  full 


Lorentz  group,  it  mist  have  the  form 


(6.56)  t  r  ,  Ia  ,  eb  ,  hc) 


where 


/  *  -  1,2,3  ) 

(  b  *  0,1,2  I 

\  c  ■  0, 1, . « . ,11  / 


1 
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■Hi 


IjL  -  t  r  £  , 

*2  *  tp  £2  > 

I,  «  tr  s? 

3 

(6.57)  •„-£•!  , 

C1  -  S.  £  S.  t 

e2-g£2€ 

b0-g-g  , 

b1  “  Ss.'S  t 

^2 

*5  -  5  •  * <*a  e)  '  \  ■  &  • 

b,  -  (£  •«)(«£  S)  ,  t6  .  (4  •  ©)(g£2£) 

b7  -  i  £(«X  £)  ,  b8  -  £  £2(«x  C)  ,  ^  -  (»X  £  ft)£2  £ 

b10  -  (i-«)  §£(gxg)  ,  bn-  (I-?)«£2(2xg) 

The  theory  presented  in  this  section  provides  a  theoretical  Justi¬ 
fication  for  the  nonrelativistic  theories  of  Toupin  [5]  and  Dixon  and 
Eringen  [6]  for  moving  electromagnetic  materials.  The  application  of 
the  second  law  of  thermodynamics  complements  the  researches  of  Jordan 
and  Eringen  [3],  [*»■],  Toupin  [5]  and  Dixon  and  Eringen  [6].  The  con¬ 
stitutive  equations  for  a  polai  able  material  (  ^  -  0  ,  J  »  0  , 

oQ  *  o)  are  equivalent  (under  the  appropriate  change  of  variables) 
to  those  of  Toupin  [5]  if  terms  of  the  order  v  /c  are  neglected. 

The  general  form  of  the  constitutive  equations  of  Dixon  and  Eringen  [6] 

are  valid  if  one  replaces  £  by  £  +  and  £  by  £  -  v-g' E 

c  c 

"Hfe  are  here  using  inks  units.  In  the  mks  system  of  units  the  ratio 

y  H 

|j)|/|  — "1  i®  tb®  eame  order  of  magnitude  as  |g|/|x  x  g|  .  The 
c 

I  x  E 

same  is  true  for  the  rations  Ifil/lx*  £l  and  lfil/1 — 5“!  • 

c 


The  specific  forms  of  the  constitutive  equations  of  Jordan  and 
Eringen  [3],  [4]  and  Dixon  and  Eringen  [6]. are  significantly  simpli¬ 
fied  by  the  second  lav.  For  example,  the  stress  tensor,  polarization 
and  magnetization  are  independent  of  the  temperature  gradient.  For 
an  isotropic  material  the  second  lav  does  not  allov  a  term  in  the 
polarization  similar  to  the  Hall  effect  in  the  current. 
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n  d* 


-D0<  V  +  +  +  %  sr  sQei  v  -  r  sr  % 

o  p 

«•»)  ,!2^ds  y°g,1r  .(Sfitfl; 


ap 


fF  .  n  **  r  „  *  n  *°i  2&  *“  •«  „ 

-  [fCL  +  no  ^  —  -  lJCa  +  no  -f*  — - 3-  * 


*P 


"7" 


>  0 


At  any  point,  x  ,  in  apace -time  the  following  quantities  can  be 
assigned  arbitrarily: 


(6.61) 


e  >  °0  -  81  »  1  ,  S 1  ,  u1  ,  d1J  ,  S1J  ,  ds1  ,  d^1  , 

DC1  ,  Da1  ,  Dd1J  (i,J  =  1,2,3) 


3y  the  argument  leading  from  (3.13)  to  (3.14 ),  a  necessary  and  suf¬ 
ficient  condition  that  the  entropy  production  (6.60)  is  non -negative 
is  that  the  following  relations  hold : ^ 


(6.62) 

f 

11 

sj0* 

(6.63) 

*0  g7  +  ^0 
*  *5 

(6.64) 

* 0 

-TT  *  n 

(6.65) 

^  -  0 

^pa 

* 


*'8  „  +  u7£  a  2 


0 


^Since  not  all  components  of  the  tensors  in  cat  (6.3d)  are  independent, 
in  deducing  (6.62)  to  (6.67)  one  should  consider  tQ  as  a  function  of 


e  ,  *r0  ,  e1  ,  s  1  ,  e  1  ,  diJ  ,  u1 
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(6.66) 

*o  sr 

a 

0  a 

(6.67) 

Ma  ’ 

-n  —  s7 

01 

(6.60) 

Plas 

91  +  JL[a  13  91  1  o 

(6.69) 

2  ..op  op 

-no  3s:  +d* 

^  ,  t  --  0 

D  D  9  d  G£ 


(6.70) 


a  * 

op  I  q  ea  .  .a 

D*  dcp  -  — 


The  following  general  statement  can  be  made  about  the  form  of  the 
constitutive  equations :  A  necessary  and  sufficient  condition  that 
the  constitutive  equations  satisfy  the  second  lav  of  thermodynamics, 
the  balance  of  moment  of  energy -momentum,  and  are  covariant  under 
proper  orthoohronous  Lorentz  trans formations  is : 
i.  The  free  energy  assumes  the  form: 


(6.71)  Vq  *  t0(e  t  i  f  ^2  9  J5^ 

vhere 

(6.72)  jx  =  £•£  ,  j2  3  &' £  •  J3  s 


ii. 

Die  entropy  ,  the  stress  tensor  t  ,  the 

polarisation  vector  *?? ,  and.  the  magnetisation 

vector  df-  are  related  to  the  free  energy  by: 

(6.73) 

\o 

1 

1 

l 

-sfe 

i 

1  ■  n 

(6*74 ) 

t  : 

2 

=  -o8  ST  £  +  i£ 

0 

(6.75) 

?  : 

dt  cfy  _ 

=  -2n  ®  -  n  °  « 

o  aJT  -  o  dJ^  ~ 

a  ,7 

(6.76) 

M 

8t  dt 

=  -an  e  -  -B 

iii. 

Die 

nonrae  chad  cal  momentum  p  is  determined  by: 

(6.77) 

£ 

K  q  -  ‘px  ??  -  Sx  ^ 

i  r. 

Die 

conduction  current  j  ,  the  heat  flov  vector 

i  ' 

and  the  dissipative  part  of  the  stress  tensor 

have  the  form: 

*  _  _  *  *  *  ^  ^  *2*  *2  *2 

(6.78)  *  -  +  k£  +  IC5«+  KkH  £  +  g  +  1^4 S  +  £  +  Kg£  £  +  ^  ? 

+  K10|  X  g  +  Kial  X  §+  K12  g  X  g  +  K15fi  XU  +  X  1  S 

*  *  *2*  ^  *2~  «.  *9„  *  *  *2* 

+  ki5§  X  4  %  +  ki63  X  4  £  +  rf^x  &  j§  +  K10^x  &‘g  +  k194  £  X  4  £ 

+  k20£  4  X  ,1  £  +  Kg^  gx  4  3  +  «22=  x  &  ^  +  k25^x  4  £ 

^  *  *  *  „  *  *  _  _  *  ^  _  *2* 

+  KgJ^X  4£  +  Kgjfi  X  4<C  +  XggfiX  4«+  Kg-^X  4  ®+Kg0#X  4  £ 


NCeMMMRfMUff  Wf«.>W 
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+  <25®  X  I2!  +  <j^x  d2!1  +  k?1£  X  l2^  +  X  42*+  <jj£ X  4*1? 

+  <<^4  £  x  4  g  +  <J5 4  ?xii+  <3^4  ?<4l+  <J74  £  *  4  £ 


+  ^  £  x  42^  ■*■  <594  £  x  I2® 


J  -  ffji-  +  ^  +  ®3©+  %4  £  +  ®j4  €  +  ag4  5+  aT4.2£  + 

+  0^^+  a10|  x  g  +  au£  x  x  8+  o^jz  x  ft  £ 


_  ,  „  „  »o*  _c 

+  al4£ X  4  £  +  a15?x  46  +  °1&e  X  4  £  +  <j17£ X  4  € 

*o  *  *  *2*  *  _  *2_  *  _  *2_, 

+  c^q-Sx  &  2?  +  0!9&  £  X  4  £  +  ct2cA  £  x  d  jg  +  a2lA  S  X  i  ^ 

+  11+  023^x  4  I  +  a24@  x  4  S  +  ff25!  xU 


*2. 


(6=79) 


*2* 


*2* 


+  o26£  xdg+  a^jff  x  4  6  +  o2gf  X  4  £  +  a29«x  4  £ 

+  oJ0®x  42C  +  a31£  x  d2€  +  ®j2£  x  I2e+  oJ5  g  x  42© 

*  *p*  *  *2*  *  <*,  *2- 

+  0344  £x  4  £  +  o}54  gx  4  £  +  ®304  5<4t 

*  *  *2„  **  *2^,  *  „  *2„ 

+  o?74  £  X  d  g  +  x  4  ®?94  §*&& 


jjfe.  -  T^g.  +  Xg4  +  Xjl2  +  \£0£+  XjgS  1+  *56®  6+  ^[^(SS+jg®  £l 
+  Xg[£«  §  +  60  jg]  +  Xjtg®  8  +  S®  g]  +  X1Q[£*(£  x  £)+(j  x  g.)0  £] 
+  ?vu[£  0(£  x  B  Mi  x®)0|]  +  X^tjg  ®(£  xS)+(gx?19  £) 

(6.8o)  +  *  £>+(£x  jg)«£]  +  \i-B?®U?x  iMSxj)  ©"fi] 

+  ®(6xSW2?xS)®«)  +  Xl6[€«(£x§M|xS)0g] 

+  \T[S0(ixg>+(£xg)®  g]  +  Xl8[|®  4£  +  4£®£] 

+  xl9I«e  ®  4  €  +  4  €1  +  ^g0[g ®  4  6  +  4  6® 6) 

+  *21tJ  ®  l2^  +  l2!  ®  J]  +  ^22^  0  42€  +  42€  ®  «] 

+  >v,j [®®42l?+  42S06i  +  *24[£®4S  +  4<Se  J) 


+ 

+ 


V,5[it®4§  +  4fi®  1]  +  *26l£8  45  +  +  ^7I)f«4l  +  ll«t] 

\tfS9 4 1  +  11®W]  v  +  4*®fi)  +  A^tl®  ¥g+  ¥g®h 


*  #o  #o  #.  „  *2  *2  „  *2*  *2* 

+  ®  4  ©+  4 €«  £l  +  4  3?+  4  fi« 41  +  *33(£®  4^  +  4^  0 €) 


,  _  *p  -fro  *o*.  *p*  *  ^  *2*  *2*  * 

+  &  ig  +  4  £<9  8]  +  *35(»®  4  ft  +  4  ft®  fil  +  ^[4  i®  4^  +  4^  x  4 

*  _  *p*  *p+  *  *  *p  *0  * 

+  *37te  §Q  O  +  O  0  4  81  +  £0  4S  +  004  a 


+  A^d  X  4  -  4  X  1]  +  X^0[|‘  X  4  -  4  X  £]  +  \-J3x  4  -  4  X  61 
+  \glix42  -  42x6]  +  \J.g  X  42  -  l2  X  C]  +  xu[»  X  42  -  42  X  «] 

+  \g[4(4  X  4  )  -  4  (e  x  4)]  +  \6(4(£x  4  )  -  4  (*x  4)] 

+  \?[4(2>x  I2)  -  42(8x  4)  ]  +  \8[I&  (4  x  4  1)  +  (1  x  4  1)  &  1) 

+  \9H& (i“x  11)  +  (ix  4  I)  0ft]  +  ^013®  »x  45)  +  (©x  4£)  »Jfi] 
+  ^(4  0  (4  x  |2|)  +  (4  x  I2!)  0  41  +  A5g[£®  (g  x  12«)  +  (£  x  42i)  ®  £3 

+  *5;[S«  (»x  4^)  +  (8  x  4%)  ®  S)  +  \*(1«  I(fc<fi)4(fcx|)  s  4] 

+  Xjjtl  0  I(6x£)4(«x4)  *  41  +  >^g[£0  4(SX0>4(£  xfi)  0  £] 

+  ^[£<8  4(€xl)4($/l)  ®£)  +  X^g[®0  4(5x4)4tex£)  ®£) 

+  ^9(8®  4(ex«)4(?/i)  «  S)  +  X6q[4  ®  (£  x  4  l)+(£  x  4  4)  0  4] 

+  A6l(4  8  (5x4  4)+(0x  44)  ®4)  +  x62[~®  (fix  4€)+(5x  4l)  ®£) 


+  A6j[£  0(1x4  S)+(4  x  4  §)  0  £1  +  ^(S  8  (4  x  4  g)+(4  xl|)«5] 
+  a65(»®  (£x  45)t(fx  4®)  ®  gi  +  yle  (g  x  l^w^x  I2!.)  ®  4) 

*  *p*  -jfp*  4-  *p  #p 

+  a67U0  (§x  4"g)+(fix  4  4)  X  4]  +  >>68(£0  (5x  4  €)+(©x  4  4)  ®£J 
+  X69[^®  (4  x  42i)+(l  x  42£)  ®  g!  +  T^S®  (4  x  4^)+(l  x  4^3)  0  & 
+  (f  x  x  t2fi)  ®  g)  +  ^(4(4  xS)«£tft  4(1  x?)l 
+  A73(4(I'x  4)  ®  fi+  5®  1(£  x  4))  +  ^(€0  (fx  I  4)+(£x  4l)®£) 
+  A^t!  @  («x  42g)+(g  x  42£)  04] 


»^gya«aw8raag=a 
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where  k^9  j  cj^,  •  ••>  9  ^ 

of  the  following  invariants) 


i  *  •  •; 


Ns  — 


functions 


(6.81) 


(6.82) 


ld  ' 
* 


II. 


Ill, 


l  > 
* 


'2 


*  * 
,  6  •  $ 


£•  &  >  §*£  ,  §■ $*  &  , 


*  *  *  ^  *  *  -  *  *  _  „  * 
&i&  ,  as  ,  %i8  ,  as  ,  £4.0  ,  US, 


U2l  ,  Sfi  ,  Sl2^  .  !?«.  ll2?  ,  f?*  , 

£4(|x  I)  ,  £ 4(6  x  £)  ,  £4(1x8)  ,  g  1(1  x *),  £4(6x8)  , 

8  4(f  x  8)  ,  £  4(1  x  Q)  ,  #4(1  x  g.)  , 

1  INg  X  1)  ,  g  ^(gx  1)  ,  1  S®(£  *  ®  ,  42(1  X  g)  ,  £  fl?(g  X  g)  , 

elhfxa) ,  £^(lx@  ,  glNlxfi) , 


,*a 


*  *P,*  *  *  *2,  ~  %  _  *2,  _  Tr 

£  4  Ce  X  4  £)  .  £4  (£  X  46.)  ,  «4  (0x  4 


me,  11),  slz(l  xk),l  ftp*  4 1)  .  slHl  X  4«)  . 
£42(gx  4  S)  ,  §42(  <£  X  4®  ,  s  FOx  4  1) 


Y»  The  heat  flux  q  ,  the  conduction  current  J  ,  and 
the  dissipative  stress  tensor  satisfy  the 
inequality: 


d 


a 

<1  ^ 


D~  -*  +^a  *  0 


which  in  particular  implies  that  when  *  0  , 


eQ  -  °  ,  =  °  ,  then 


l 


10$ 


(6.85)  Dt°*  -  0  ,  <ia  -  0  ,  f  -  0 

Th ese  ]  vaults  follow  directly  from  equations  (6.62)  to  (6. 70)  and 
invariance  requirements.  In  particular ,  (6.64)  and  (6.65)  imply  that 

+o  *  +o(e  '  n0  »  Sa  >^a  »  V 

and  in  order  that  be  an  invariant  under  the  proper  orghochronous 
Lorentz  group,  it  must  have  the  functional  form  (6.71).  Equations 
(6.73)  to  (6.76)  follow  directly  from  (6.62),  (6.69),  (6.66),  (6.67) 
and  (6.71).  Equations  (6.65)  and  (6.68)  are  identically  satisfied 
(cf.  (6.71),  (6.75)  and  (6.76)).  Ike  lengthy  expressions  (6.78)  to 
(6.60)  are  a  straightforward,  application  of  the  remarks  in  the  appendix. 

Ike  foregoing  theory  of  viscous  fluids  is  so  complicated  that  in 
the  generality  presented  in  this  section  it  has  little  practical  use. 
However,  it  provides  a  foundation  for  several  approximate  theories. 

In  particular,  the  class  of  perfect  fluids  is  a  special  case  of 
this  theory. 


Perfect  Fluids 

For  a  perfect  fluid,  the  entropy  production  vanishes.  That  is: 


(6.84)  q.0>Dt-O,j.O 


Hint  is,  for  a  perfect  fluid,  the  saterials  are  described  by: 

i.  The  free  energy  has  the  form  (6.71)* 
ii.  The  entropy  ^  ,  the  polarization  vector  f?, 
and  the  magnetization  vector  are  given  by 
.  (6.73),  (6.75)  and  (6.76)  respectively  and  the 

stress  tensor  is  determined  from: 

2 

(6.85)  i.  «  -  nQ  j— •  & 

ill.  Hae  noon*  chani  cal  nomentum  p  is  determined  by: 

(6.86)  p  *  -  §  -  j§x  t!? 

The  theory  of  perfect  electromagnetic  fluids  presented  in  this 
section  ia  comparable  to  the  theory  of  Penfield  and  Eaus  [36].  That 
the  polarization  and  magnetization  are  derived  from  the  free  energy  is 
a  consequence  of  the  second  lav  of  thermodynamics  in  the  theory  of 
this  section.  In  the  theory  of  Penfield  and  Baus  the  polarization 
and  magnetization  are  deduced  from  a  variational  principle.  The  major 
difference  between  the  work  of  this  article  and  that  of  Penfield 
and  Ha us  is  i r  vne  difference  in  the  Amperian  model  and  the  magnetic 
dipole  model  for  the  magnetic  term  in  the  interaction  of  electro¬ 
magnetic  fields  on  matter. 


Concluding  Re  Mirks 


Hie  relativistic  theory  of  electromagnetic  materials  presented  in 
this  paper  provides  an  insight  into  the  continuum  behavior  of 
electromagnetic  interactions  with  matter.  The  coupling  of  mechanical, 
electromagnetic  and  thermal  phenomena  into  one  ‘theory  has  resulted  in 
e,  complicated  system  of  partial  differential  equations.  The  complete 
physical  and  mathematical  nature  of  these  equations  in  the  generality 
treated  in  the  previous  sections  is  beyond  the  scope  of  our  present 
knowledge,  Hie  nonlinearity  of  the  equations  makes  any  solution  to 
boundary -value  problems  difficult. 

Considerable  insight  into  the  nature  of  a  system  of  partial  dif¬ 
ferential  equations  can  be  obtained  from  an  investigation  of  the  propa¬ 
gation  of  waves  and  singular  surfaces.  This  seems  to  be  a  fruitful 
class  of  problems  which  can  be  treated  relativistically.  Some  progress 
in  this  direction  has  already  been  made  and  will  appear  in  %  later 
work. 

T1  "nigh  the  present  theory  is  capable  of  describing  many  physical 
phenomena  (polarization,  magnetization,  heat  conduction,  piezoelectricity, 
thermoelectric  effects,  to  name  a  few)  such  physical  phenomena  as 
viscoelasticity,  optical  activity  and  gyrotropic  effects  are  definitely 
excluded.  Simple  theories  of  viscoelasticity  (the  Kelvin-Voigt, 

Maxwell  and  other  higher  order  rate  theories )  can  be  formulated  by 
an  appropriate  change  in  constitutive  equations.  Hie  treatment  of 
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: 


heredity  requires  the  more  difficult  study  of  functional  constitutive 
equations.  A  description  of  optical  activity  and  gyrotrqpic  effects 
will  probably  entail  a  reformulation  of  the  balance  equations  to 
include  spin,  couple  stresses,  and  quadrupoles .  Whether  ferromagnets 
and  electrets  are  described  by  the  theory  presented  in  this  article 
cannot  be  answered  until  a  deeper  physical  study  of  these  phenomena 
is  undertaken. 

Finally  it  must  be  mentioned  that  for  accelerating  frames 
and  curved  spaces  where  the  special  theory  fails  new  unified  theories 
are  needed  employing  the  fundamental  ideas  of  the  general  theory 
of  relativity.  Such  a  grandiose  plan  presently  is  out  of  our  reach. 


110 


BIBLIOGRAPHY 


1.  Toupin,  R.  A.,  "The  Elastic  Dielectric,"  J.  Rat.  Blech.  and. 

Anal.  2,  849-915  (1956). 

2.  Eringen,  A.  C. ,  "On  the  Founds ti oni  of  Electroelastostatics, N 
Int.  J.  &i£.  Sci.,  X,  127-154  (1965). 

3.  Jordan,  H.  F.  and  Eringen,  A.  C.,  "On  the  8tatic  Vonlinear 
Theory  of  Elec  tr^agne  tic  Theraoelastic  Solids  -  I,"  Int.  J. 

Enft.  Sci.,  g,  36-96  (1964). 

4.  Jordan,  If.  F.  and  Eringen,  A.  C. ,  "On  the  Static  nonlinear 
Theory  of  Electromagnetic  Thermoelastic  Solids  -  II,”  Int. 

J.  Bn*.  Sci.,  g ,  97-114  (1964).  “ 

5.  Tcwpin,  R.  A.,  "A  Dynamical  Theory  of  Elastic  Dielectrics,” 

Int.  J.  En^.  Sci.,  i,  101-126  (1963). 

6.  Dixon,  R.  C.  and  Eringen,  A.  C.,  "A  Dynamical  Theory  of  Polar 

Elastic  Dielectrics  I  and  II,"  Int.  J.  Ifcg.  8ci.,  359-598 

(1965). 

7.  Eringen,  A.  C.,  lonlinear  Theory  of  Continuous  Media,  McGraw- 
Hill,  lev  York,  19551 

8.  Truesdell,  C.  and  Toupin,  R.  A.,  "The  Classical  Field  Theories," 
Handbuch  der  Physlk,  III/l,  8pringer-Verlag,  Berlin,  i960. 

9.  Eringen,  A.  C.,  "A  Unified  Theory  of  Thermomechanical  Material," 
Int.  J.  Eng.  Sci.,  4,  (1966). 

10.  Bressan,  A,,  "Cinesatlca  del  Sistemi  Ccntinui  in  Relativista 
Gene  rale,"  Annali  di  Meicnatica  Pure  ed  Appllcata,  62.,  Series  4, 
99-lUB  (196577 

11.  Toupin,  R.  A.,  "World  Invariant  Kinematics,"  Arch.  Rat.  Mech. 
Ami.,  I,  181-211  (1957). 

12.  Synge,  J.  L.,  Relativity:  The  Special  Theory  1  forth -Holland 
Publishing  Co.,  Amsterdam,  1965* 

13.  Lichnerovicz,  A.,  Theories  Relatives tes  de  la  Gravitation  et 
de  1  'Elccti omagpdtlsae ,  Masson  et  Cle,  Paris,  1955. 

14.  Miller,  C.,  The  Theory  of  Relativity,  Oxford,  London,  i960. 


Ill 


15. 

16. 

17. 

18. 

19. 

20. 

21. 

22. 

23. 

24. 

25. 

26. 

27. 

28. 
29. 


Landau,  L.  D.  and  Llfshitz,  E.  M. ,  Fluid  Mechanic* ,  Fsrgamon 
Press,  London,  1999.  -  - 

Van  Dantzig,  D. ,  "On  Relativistic  Thermodynamics,”  Kon.  Red. 

Acad.  Van  Wetenschappen  Proc.,  42,  601-607  (1939). 

Bckart,  C.,  "The  Thermodynamics  of  Irreversible  Processes  III. 
Relativistic  Theory  of  Sisple  Fluid,"  Phys.  Rev.,  $&,  919-924  (1940). 

Bergmann,  P.  G. ,  Introduction  to  the  Theory  of  Relativity, 

Prentice -Ball,  Englewood  Cliffs,  f.  T. ,  1942. 


Thomas,  T.  T.,  Concepts  from  Tensor  Analysis  and  Differential 
Geometry  (2nd  Edition; 7  Academic  ^re os,  I.  Y. ,  1963. 


Edelen,  D.  G.  B. ,  The  Structure  of  Field  Space,  University  of 
California  Press,  Berkeley  and  Los  Angeles,  1962. 


Bogoliubov,  H.  If.  and  Shirkov,  3).  V. ,  Introduction  to  the  Theory 
of  Quantized  Fields,  Inter  science.  If.  T.,  1959- 

Papapetrou,  A.,  "Jfon-Syametric  Stress -Energy-Momentum  Tensor  and 
Spin  Density,”  fhll,  Mag.,  40,  937-946  (1949). 

Bragg,  L.  E.,  "On  Rslatlvietic  Worldlines  and  Motions  and  on 
Ron-Sentient  Response,"  Arch.  Rat.  Mech.  Anal.,  l8,  127-166  (1963). 

Bressan,  A.,  " The rmodi name ia  e  Magneto- Vis co -Blast iciti  con 
Deformazioni  Finite  in  Re  la ti vita  Gene rale,"  Rend,  del  8em. 

Mat,  della  Univ.  di  Palova,  J4,  1-73  (1964). 

Pham  Mau  Quan,  "Sur  une  Theorie  Relativists  des  Fluides  Therrno- 
dynamiques , "  Annali  di  Matematica  Pura  ed  Applicata,  38,  Series  4, 
121-203  (1955T 

Fano,  R.  K. ,  Chu,  L.  J.  and  Adler,  R.  B.,  Electromagnetic  Fields, 
Energy,  and  Forces ,  John  Wiley  and  Sons,  Inc.,  I.  l7,  I960. 


To  1 man,  R.  C.,  Relativity,  Thcj 
at  the  Clarendon  Press,  1934. 


.ca  and  Cosmology,  Oxford 


Landau,  L.  D.  and  Lifshitz,  E,  M.,  The  Classical  Theory  of  Fields 
Pergamon  Press,  London,  19 62.  “  ^  ^  ~ 


Eisenhart,  L.  P.,  Continuous  Group  of  Transformations ,  Princeton 
University  Press,  Princeton,  H..J.7  I9S37"lbover  5Pukl.,  I.Y.,  1963). 


112 


30.  Weyl,  H. ,  Lpace -Time -totter t  Dover  Publications,  I.  T.,  1982. 

31.  d'Abro,  A.,  The  Evolution  of  Scientific  Thought  from  Hevton  to 
Einstein,  Dover  Publi cations,^!.  ¥., 

32.  Noll,  W.,  "The  Foundations  of  Classical  Mechanics  in  the  Light 
of  Recent  Advances  in  Continuum  Mechanics,"  in  The  Axiomatic 
Method  with  Special  Reference  to  Geometry  and  Physics,  edited 
by  Henkin,  L. ,  Suppes,  P.,  andTJar^ki,  A"; *l?orth -Holland 
Publishing  C/onqjany,  Amsterdam,  1939* 

33-  Smith,  F.  G.,  "On  Isot;  opic  Integrity  Bases,"  Arch.  Rat.  Mach. 
Anal.,  l8,  282  (1965) 

34.  Post,  E.  J. ,  Formal  Structure  of  Electromagnetics,  Bor th -Holland, 
Amsterdam,  19&2. 

33*  Lorentz,  H.  A.,  The  Theory  of  Electrons,  Dover  Publications,  I.  T. , 
1952. 

36.  Penfield,  P.  and  Baus,  H.  A.,  "Energy -Momentum  Tensor  for  Electro¬ 
magnetic  Fluids, "  in  Recent  Advances  in  Engineering  Science , 
edited  by  Eringen,  A.  C.,  Gordon  and  Breach  Science  Publishers, 

K.  y.,1966. 

37.  Wine  man,  A.  S.  and  Pipkin,  A.  C.,  "Material  Symmetry  Restrictions 
on  Constitutive  Equations,"  Arch.  Rat.  Mech.  Anal.,  17,  l84  (1964). 

3fl.  Truesdell,  C.,  "General  and  Exact  Theory  of  Waves  in  Finite 
Elastic  Strain,"  Arch.  Rat.  Mech.  Anal.,  8,  263  (196l). 

39.  Weyl,  H. ,  The  Classical  Groups ,  Their  Invariants  and  Representa¬ 
tions,  Princeton  University  Press,  Princeton,  ft.  ,  194b. 

40.  Spercer,  A»  J.  M.  and  Rivlin,  R.  S.,  "The  Theory  of  Matrix 
Polynomials  and  its  Application  to  the  Mechanics  of  Isotropic 
Continua,"  Arch.  Rat.  Mech.  Anal.,  2,  309  (1959). 

41.  Spencer,  A.  J.  M,  and  Rivlin,  h.  8.,  "Isotropic  Integrity  Bases 
for  Vectors  and  Second  Order  Tensors,"  Arch.  Rat.  Mech.  Anal. ,  9, 

43  (1962).  ~~  ™ 

42.  Spencer,  A.  J.  M.,  "Isotropic  Integrity  Bases  for  Vectors  and 
Second  Order  Tensors,"  Arch.  Rat.  Mech.  Anal.,  18,  51  (196?). 


1X5 


APPENDIX  OH  POLYNOMIAL  INVARIANTS  OP  VECTORS 
AND  TEKSGRS  POR  THE  LGHENTZ  GROUP 


H -Vectors  VQ^  (i  -  1, *♦*,!) 

Tbit  polynomial  invariants  for  the  Loren  tz  group  can  be  deduced 
from  the  polynomial  invariants  of  the  four  dimensional  orthogonal, 
group  (see  Weyl  [59,  p.65]).  If  f  is  an  even  invariant  polynomial 
of  N-vectors  (i  -  1, »..,N)  ,  then  f  is  a  polynomial  in  the 

scalar  products 


(A.l)  '/*  Va(i)  Vp(j)  (i  <  J  ;  i,4  -  i,...,*) 

If  f  is  an  odd  invariant  polynomial  of  E-vectors  ,  f  is 

a  sum  of  terms  of  the  form: 


(A.2) 


g 


where  g  is  an  even  polynomial  and 

ix  <  i8  <  ij  <  i4  (ix,  i2»  i4  -  X, •••,») 

The  following  special  case  is  often  encountered  in  physical  situations. 


11* 


(J ) 

N-Vectors  V  1  J  and  a  Tine -Like  Unit  Vector  ug 
(u^  ua  "  -l)(i  » 

In  many  physical  situations  the  vorla  velocity  vector  ua  is 
included  in  the  group  of  vectors  considered  in  Case  1.  The  world 
velocity  is  time -like 


(A.3)  ua  ua  -  -1 


Any  polynoadal  in  V  ^  and  ua  is  a  polynomial  in  the  vectors 


5  U) 

a 


and  uQ  and  the  scalar  products 


(A.*)  Y(l)  •  V  (1)  u“  (i-1, 


*  (i) 

where  V^'  '  is  defined  as 


(A.5)  va(l)  •  sj*  Vp(i)  i  S £ 


bci  +uau^  '  V  ^  ua 
a  a  'a 


This  follows  from  the  decomposition 


(A. 6)  Ya(i)  -  VQ(i)  -  v(i)  ua 


Since  the  v^  are  scalar  invariants,  one  need  only  consider  a 


polynomial  invariant  in 
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(A. 7) 


$  <*> 

a 


u 


a 


By  applying  the  results  of  Case  1  to  the  set  (A. 7)  and  observing 
(A. 8)  u“va(i)  =  0  ;  uaua  .  -1 


aid. 


(A. 9) 


€<W  J  (ii)  $  (12)  J  U5)  J  Uu) 

a  p  7  b 


0 


#  /  i  } 

one  obtains:  An  even  invariant  polynomial  f  in  the  vectors  V  v  '  , 
ua  is  a  polynomial  in  the  invariants 


(A. 10) 


J*  1  (i)  $  (J) 

’  a  vp 


(1S3HJ-  1, 


and  an  odd  invariant  polynomial  f  is  a  sum  of  terras  of  the  form: 

(a. ii)  8**$oUi>  y*a)  ^5>  u6  8 

(lx  <  i2  <  lj  ;  i1(  i2,  ij  -  1....H) 


where  g  is  an  even  invariant  polynomial. 

Ihe  identity  (A. 9)  follows  from  the  observation  that  from  (A. 8) 
the  fourth  row  of  the  determinant  on  the  left  hand  aide  of  (A. 9)  is 
a  linear  coabi2irti.cn  of  the  first  three  rows. 


The  procedure  of  decomposing  tensors  into  spatial  tensors  and 
scalars  will  be  used  in  the  next  section.  It  is  the  key  to  using 
the  results  of  the  three  dimensional  orthogonal  group  for  the  integrity 
bases  of  the  Lorentz  group. 


M-Syraoetrlc  Tensors  ,  If -Vectors  VQ^ 

and  a  Tine -Like  Unit  Vector  (u^u®  -  -1) 

Any  polynomial  in  ,  and  34  is  a  polynomial  in  the 


tensors 

|(R) 

,  the  vectors 

*(i) 

2  ,  h  ,  14  and  the  scalars 

(A.12) 

t(r) 

■  >£’ •? 

(R  »  1,...,M) 

(A  15) 

VU) 

■  \w »° 

(i  -  1,... 

where 


(A.14)  TgJ  «  -  rg> 


(A.15)  Va(i)  ■  sj3  vRCi) 


(A. 16)  Aa(R)  *  S^7  T^R)  U5 
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This  follows  froB  (A. 6)  and  the  identity 


/*  I'M  *(R)  2i<R)  *  (*)  *  (*)  mW 

(A.17)  -  Aav  -  Apv  ua  ♦  Tl  n 


au0 


It  should  he  noted  that 


(A.l8)  *  0 

up 


(A.19)  Aa(R)  u°  -  0 


(A. 20)  Vff^  u°  *  0 


Thus  one  can  consider  a  polynomial  in  the  tensors  and  vectors 


(A.21)  T^}  ,  vjk)  ,  ua  (R  .  1,...,M  ;  k  »  1,...,M+  I) 


idle  re 


(A, 22) 


*(IHM)  i#  A  (R) 
a  a 


with  the  coefficients  polynomials  in  and  . 

The  coefficients  of  an  invariant  polynomial  must  he  isotropic 

tensors  for  the  Lorentz  group.  This  implies  that  they  are  sums  of 

as  aSr 5 

terms  which  are  products  of  y  and  c  By  noting  that 


U8 


(A. 23) 

<v* 

“0  ' 

o  , 

<V*  ua 

-  0 

(A.a*) 

Va(i)ua  .  ( 

) 

(A. 25) 

‘  0076 

i%r 

$  (i) 
p 

(‘v*  VJ) 

<V  V1 

(A. 26) 

€  0076 

V0(k) 

*  (t) 

.  0 

(A. 27) 

*  0076 

(V* 

-  0 

an  invariant  polynomial  in  J  >  )L  and  must  be  a  polynomial  in 


(A. 38)  V,<‘>  (i,)»  v  <J>  ,  t» 


«•»>  < «,  (V®  V°  <V  V"  'a1'1 


(A. JO)  •«  tf/*  V‘> 


*  *  *  .  *  * 

where  %  >  2^j  >  %  ,  Zp  >  8104  2g  fcre  nBtrix  product*  in 

*(p)  *  *(r} 

X  and  x  *  are  matrix  product#  in  J  '  .  The  proof  of  (A. 23) 
and  (A. 24)  follows  from  the  definition#  of  and  ^  .  Identities 

(A. 25),  (A. 26)  and  (A. 27)  are  deduced  by  observing  that  the  fourth  row 
of  the  determinant#  i#  a  linear  combination  of  the  first  three  (cf.  (A.l8) 
to  (A. 20)  ).  Here  c  ^ 


is  defined  by: 


•  I 


i 

t 
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^A'51^  *  C&7  *  *0076  U 


The  following  identities  are  useful  to  reduce 

(A. 28) 

to  (A. 30) 

Vs 

Vs 

Vs 

Vs 

Vs 

Vs 

Vs 

Vs 

(A,52)  \Wi 

Vs7S6S  "  - 

Vs 

Vs 

Vs 

Vs 

Vs 

Vs 

Vs 

Vs 

From  (A. 32)  one  can 

show  that 

V 

P2 

'\ 

(A. 33)  «  1  1  1  1 

Vs7S&S  “  - 

*\ 

A 

7S 

*\ 

61  &2 

By  dotting  (A. 32)  with  u  u  and  adding  it  tines  the  fourth  row 


to  the  first  row,  uft  tines  the  fourth  row  to  the  oecond  row,  u 

pi  7i 

tines  the  fourth  row  to  the  third  row,  one  obtains 
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By  dotting  the  identity 
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with  u  u  and  multiplying  the  fourth  row  hy  u  and  adding  it 

1 


to  the  first  row,  etc.;  one  can  deduce  that 
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(A. 36) 
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From  (A.  3*0  it  follow*  that 
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If  one  define*  the  skew -symmetric  tensor  by 


(A.  38) 
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(A. 39) 


v<i>uP  -  0 


From  (A. 38)  it  follow*  that 

“•"«  Vi!  -  s  !(1)67 


By  substituting  (A.40)  into  (A. 28)  to  (A. 30)  it  is  seen  that  tbs 
polynomial  invariants  of  N  symmetric  tensor,  I-vectors  and  a 
time -like  unit  rector  reduces  to  the  set  of  invariants  (A.  12)  and 


(A.  13)  and  the  invariants  of  N  syne  trie  and  N  +  V  skev-synttric 
space -like  tensors  which  are  polynomials  of  degree  three  or  less  in 
the  skev-symeetric  tensors.  This  is  the  starting  point  of  the  work 
of  Rivlin,  Spencer  and  Snith  (cf.  [40],  [4l],  [42],  [33])  for  three 
dimensional  matrices. 

From  (A. 36)  space -like  tensors  satisfy  an  equation  vhlch  in 

aatrix  notation  is  identical  to  the  three  dimensional  Clayley -Hamilton 

*^1^2 

theorem.  If  one  multiplies  (A.  36)  by  a  a  a  and 
expands  the  determinant  the  following  result  is  obtained: 


(A.lfl) 


ft3  -  <tr  |  t(tr  If  .  tr  £2]  J  -fc  [tr  V 

+  2tr  ft5  -  3(tr  ft2)(tr  ft)]  £  •  0 


*03 

where  a  satisfies 


(A. 42)  EPfa  Pb  .  2* 


Since  it  is  only  this  theorem  and  several  properties  of  the  trace  of 
products  of  matrices  that  are  employed  by  Rivlin,  Spencer  and  Smith 
in  their  series  of  papers  on  the  reduction  of  polynomials  in  3  x  3 
matrices,  their  results  (cf.  [40],  [4i],  [42],  [33])  bold  for  space 


r 


1*3 


i 


Ilia  tensors  end  space-like  vector*  under  the  Lorentz  group.  (In 
particular  all  relations  listed  by  Spencer  [42,  p.54J  are  valid  tor 
•pace -like  tensers.  Only  these  relations  are  ©ployed  by  these  authors 
for  the  reduction  of  polynomial  Invariant  of  3  x  3  matrices  under  the 
orthogonal  group.) 

The  integrity  bases  listed  by  Spencer  [42]  and  Smith  (33 3  are 
minimal  in  the  sense  that  no  invariant  listed  by  them  is  a  polynomial 
in  the  invariants  in  their  list.  The  integrity  bases  for  space -like 
tensor  derived  from  their  results  is  also  minimal  in  the  above 
sense.  If  it  were  not,  it  would  imply  that  their  result  was  not 
minimal.  Smith  has  proven  that  their  results  (cf.  [36])  are  minimal . 

The  results  of  Vlceman  and  Pipkin  [37]  for  continuous  tensor 
functions  of  tensors  for  the  orthogonal  group  In  three  dimensions  hold 
for  a  space -like  tensor  function  of  uQ  and  four -tensors  (even  though 
the  Lorentz  group  is  not  compact).  The  fact  that  a  polynomial  basis 
in  ua  and  four  tensors  is  a  functional  basis  can  be  proven  by  their 
method  or  by  observing  that  if  it  were  not  a  contradiction  would 
arise  with  their  results.  By  transforming  to  a  frame  in  which 
ua  *  observing  that  the  resulting  tensor  must  be 

invariant  under  the  three  dimensional  orthogonal  group  one  can  see 
that  the  above  remaik  as  to  the  form  of  a  continuous  invariant  function 
must  hold  for  space -like  tensors. 
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